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Abstract. In this paper we study the metastable behavior of one of the simplest 
disordered spin system, the random field Curie-Weiss model. We will show how 
the potential theoretic approach can be used to prove sharp estimates on capaci- 
ties and metastable exit times also in the case when the distribution of the random 
field is continuous. Previous work was restricted to the case when the random 
field takes only finitely many values, which allowed the reduction to a finite di- 
mensional problem using lumping techniques. Here we produce the first genuine 
sharp estimates in a context where entropy is important. 
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1. Introduction and main results 



The simplest example of disordered mean field models is the random field Curie- 
Weiss model. Here the state space is Sn = {— 1, 1}^, where iV is the number of 
Ph ; particles of the system. Its Hamiltonian is 

> . AT I 1 \ 

|: H N [u)(tr) = -- (l^X>) "E*. (l.D 

where A = {1, . . . , iV} and hi, i 6 A, are i.i.d. random variables on some probability 
space (fi,^ 7 , P/i). For sake of convenience, we will assume throughout this paper 
that the common distribution of h has bounded support. 

The dynamics of this model has been studied before: dai Pra and den Hollan- 
der studied the short-time dynamics using large deviation results and obtained the 
\q '. analog of the McKeane-Vlasov equations [fl6ll . Mathieu and Picco Ifl5ll and Fontes, 

Mathieu, and Picco [fl2ll . considered convergence to equilibrium in a particularly 
simple case where the random field takes only the two values ±e. Finally, Bovier 
et al. ^ analyzed this model in the case when h takes finitely many values, as 
an example of the use of the potential theoretic approach to metastability. In this 
article we extend this analysis to the case of random fields with continuous distri- 
butions, while at the same time improving the results by giving sharp estimates of 
transition times between metastable states. 

The present paper should be seen, beyond the interest presented by the model as 
such, as a first case study in the attempt to derive precise asymptotics of metastable 
characteristics in kinetic Ising models in situations where neither the temperature 
tends to zero nor an exact reduction to low-dimensional models is possible. While 
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the RFCW model is certainly one of the simplest examples of this class, we feel that 
the general methodology developed here will be useful in a much wider class of 
systems. 

1.1. Gibbs measure and order parameter. The static picture. The equilibrium 
statistical mechanics of the RFCW model was analyzed in detail in 0]] and lfl3H . 
We give a very brief review of some key features that will be useful later. As usual, 
we define the Gibbs measure of the model as the random probability measure 

2 -N e -l3H N [u](a) 

/49,ivM(<7) = — y — n — , (1-2) 
where the partition function is defined as 

Zf,, N [u>] = E ff e"^ MW = 2~ N e"^ H{,) . Q-3) 



We define the total magnetization as 



(-)-4$> (L4) 



m N ^ - N 

i6A 



The magnetization will be the order parameter of the model, and we define its 
distribution under the Gibbs measures as the induced measure, 

Qp,N = (J>p,N ° «V\ (1-5) 

on the set of possible values = {—1, —1 + 2/N, . . . , 1}. 
Let us begin by writing 

Z/3,N[u]Qi3 tN [u](m) = exp ^-^- m ^ ^/3,ivM( m ) 0-6) 

where 

\ ieA / 

For simplicity we will in the sequel identify functions defined on the discrete set 
T N with functions defined on [—1, 1] by setting /(m) = f([2Nm]/2N). Then, for 
m G (—1,1), Zjj(m) can be expressed, using sharp large deviation estimates [10], 
as 

Zl N [u]{m) = 7 Ni U U (1 + o(l)) , (1.8) 

where o(l) goes to zero as N f oo. This means that we can express the right-hand 
side in (II. 6D as 



^MQ A ivM(m) = Y^«^exp(-^/3F AJ v[a;](m)) (l+o(l)), (1.9) 

where 

F Aiv [a;](m) ee -^m 2 + ^[^(m). (1.10) 
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Here I N [u](y) is the Legendre-Fenchel transform of the log-moment generating 
function 



U N [u](t) =. y 



^lnE^exp^p^ (1.11) 
— In cosh (t + Phi) . 



N . , 

Above we have indicated the random nature of all functions that appear by making 
their dependence on the random parameter uj explicit. To simplify notation, in the 
sequel this dependence will mostly be dropped. 

We are interested in the behavior of this function near critical points of F^ N . An 
important consequence of Equations (II .61) through (II. IIP is that if m* is a critical 
point of Fp )N , then for \v\ < N~ 1/2+s , 

g ft „(m-) = exp ("T a{m )v ) {1+ o(1)) ' (L12) 

with 

a(m*) = Fp N {m*) = -1 + (3~ l I^(m*). (1.13) 
Now, if m* is a critical point of F^jv , then 

m * = (3~ l I' N {m*) = p'H*, (1.14) 

or 

p m * = I' N (m*) = t*. (1.15) 
Since I N is the Legendre-Fenchel transform of Un, I' n (x) = U'^ l (x), so that 

m* = U' N (Pm*) = 1 ^tanh(/3(m* + hi))). (1.16) 

Finally, using that at a critical point, I'^im*) = u „ l u*y we get the alternative 
expression 

GK) = ~ 1 + PUUPm*) =~ 1+ A EieA (i_ tanh ^ (m * + ^))- (L17) 

We see that, by the law of large numbers, the set of critical points converges, P^- 
almost surely, to the set of solutions of the equation 

m* = E h tanh(/3(m* + h)) , (1.18) 

and the second derivative of Fp )N (m*) converges to 

lim F"(m*) = -1 + -. \ r . (1.19) 

^v-oo P> NK ' pE h (1 - tanh 2 (/3(m* + h))) 

Thus, m* is a local minimum if 

PE h (1 - tanh 2 (/3(m* + h))) < 1, (1.20) 

and a local maximum if 

pE h (1 - tanh 2 (/3(m* + h))) > 1. (1.21) 

(The cases where pE h (l — tanh 2 (/3(m* + h))) = 1 correspond to second order 
phase transitions and will not be considered here) . 
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Proposition 1.1. Let m* be a critical point of Qp^N- Then, F h -almost surely, for all 
but finitely many values of N, 

7 n ( *\ exp(-PNFp tN (™*)) (l + o(l)) n 99 > 

Z(3,NQp,N{m ) = (1.22) 

|E (1 - tanh 2 (/?(m* + h))) \ 

with 

Ff3, N (m*) = ~ J^/Z lncosh + ^)) ■ Cl-23) 

From this discussion we get a very precise picture of the distribution of the order 
parameter. 

1.2. Glauber dynamics. We will consider for definiteness discrete time Glauber 
dynamics with Metropolis transition probabilities 

Pn[u]((t,v') = ^exp(-^[H N [cu](a')-H N [u](a)] + ), (1.24) 
if a and a' differ on a single coordinate, 

p N [cu]{a,a) = 1 - ^ ex P (-/9[fljvM(0 - #ivM0x)]+) , (1.25) 

and pat(ct, cr') = in all other cases. We will denote the Markov chain correspond- 
ing to these transition probabilities a(t) and write F u [u] = P„, for the law of this 
chain with initial distribution v, and we will set P CT = P <5ct . As is well known, this 
chain is ergodic and reversible with respect to the Gibbs measure pp^ N [u], for each 
u. Note that we might also study chains with different transition probabilities that 
are reversible with respect to the same measures. Details of our results will depend 
on this choice. The transition matrix associated with these transition probabilities 
will be called P N , and we will denote by L N = P N — 1 the (discrete) generator of 
the chain. 

Our main result will be sharp estimates for mean hitting times between minima 
of the function Fp^im) defined in (11.10D . 

More precisely, for any subset A c Sn, we define the stopping time 

t a = infO > 0\a(t) E A}. (1.26) 

We also need to define, for any two subsets A,Bc Sn, the probability measure on 
A given by 

va,b{?) = v=; / vip r — Z — r (1-27) 

E CTe A^(o-)P CT [r B < T A ] 

We will be mainly concerned with sets of configurations with given magnetization. 
For any / e rv, we thus introduce the notation S[I] = {a e Sn ■ rn N (a) G 1} and 
state the following: 

Theorem 1.2. Assume that (3 and the distribution of the magnetic field are such that 
there exist more than one local minimum of Fp iN - Let m* be a local minimum of Fp^ N , 
M = M(m*) be the set of minima of Fp )N sucn that Fp^{m) < Fj3 t N(m*), and z* be 
the minimax between m and M, i.e. the lower of the highest maxima separating m 
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from M to the left respectively right. Then, F h -almost surely, for all but finitely many 
values of N, 

^smu^siM] = exp(PN[Fw{J)-Fw(m?)]) (1.28) 

2nN I (3E h (l-ttmh 2 (l3(z* + h)))-l 

X- 7\ 7 g r- (1 + 0(1)) , 

(3\li\ y 1 - /3E h (1 - tanh 2 {(3{m* + h))) 
where 7 X is the unique negative solution of the equation 

- tanh(/3(z* + h))) exp (-2/3 [z* + h] 



exp (-20[z*+h] + ) ^ 



1. (1.29) 



/3(l+tanh(/3(z*+/i))) 

Note that we have the explicit representation for the random quantity 

F p , N (z*) - F PfN {m*) = ( ^* )2 ~ (m * )2 (1.30) 



~ m 2 [ lncosh + ^)) - lncosh (/3(m * + 



The proof of this result on mean transition times relies on the following result 
on capacities (for a definition see Eq. (12.51) in Section 2 below). 

Theorem 1.3. With the same notation as in Theorem \l.2\ we have that 

7 (Qt *1 cr^n IN exp (-PNF PiN (z*)) (l+o(l)) 
Z^jvcap J, 5[MJ) = — / = = (1.31) 



/3Eft (1 - tanh 2 ((3(z* + h))) - 1 



The proof of Theorem 11.31 is the core of the present paper. As usual, the proof 
of an upper bound of the form (I1.31D will be relatively easy. The main difficulty is 
to prove a corresponding lower bound. The main contribution of this paper is to 
provide a method to prove such a lower bound in a situation where the entropy of 
paths cannot be neglected. 

Before discussing the methods of proof of these results, it will be interesting to 
compare this theorem with the prediction of the simplest uncontrolled approxima- 
tion. 

The naive approximation. A widespread heuristic picture for metastable behavior 
of systems like the RFCW model is based on replacing the full Markov chain on S N 
by an effective Markov chain on the order parameter, i.e. by a nearest neighbor 
random walk on r w with transition probabilities that are reversible with respect 
to the induced measure, Qp^. The ensuing model can be solved exactly. In the 
absence of a random magnetic field, this replacement is justified since the image 
of a(t), m(t) = mjv(c(t)), is a Markov chain reversible w.r.t. Qp^', unfortunately, 
this fact relies on the perfect permutation symmetry of the Hamiltonian of the 
Curie-Weiss model and fails to hold in the presence of random field. 
A natural choice for the transition rates of the heuristic dynamics is 

r N [u)(m,m') = — 1 V V p N [uj)(a,a'), (1.32) 

a:m^(a)=m o" :mjv(o" )=m 

which are different from zero only if m' = m ± 2/N or if m = m'. The ensuing 
Markov process is a one-dimensional nearest neighbor random walk for which most 
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quantities of interest can be computed quite explicitly by elementary means (see 
e.g. ITT71 I3T). In particular, it is easy to show that for this dynamics, 

E v s[m *),siM] T siM] = exp(pN[F 0tN (z*)-F fi)N (m*)]) 



2ttN (3E h (l-tznh 2 (f3(z* + h)))-l 
X [3\a(z*)\\j l-(3E h (l-tanh 2 ((3(m* + h))) [ ° Uj ' 

where a(z*) is defined in Q1.19D . 

The prediction of the naive approximation is slightly different from the exact 
answer, albeit only by a wrong prefactor. One may of course consider this as a 
striking confirmation of the quality of the naive approximation; from a different 
angle, this shows that a true understanding of the details of the dynamics is only 
reached when the prefactors of the exponential rates are known (see |fT4ll for a 
discussion of this point) . 

The picture above is in some sense generic for a much wider class of metastable 
systems: on a heuristic level, one wants to think of the dynamics on metastable 
time scales to be well described by a diffusion in a double (or multi) well poten- 
tial. While this cannot be made rigorous, it should be possible to find a family of 
mesoscopic variables with corresponding (discrete) diffusion dynamics that asymp- 
totically reproduce the metastable behavior of the true dynamics. The main mes- 
sage of this paper is that such a picture can be made rigorous within the potential 
theoretic approach. 

Acknowledgments. The authors thank Alexandre Gaudilliere, Frank den Hollan- 
der, and Cristian Spitoni for useful discussions on metastability. 

2. Some basic concepts from potential theory 

Our approach to the analysis of the dynamics introduced above will be based on 
the ideas developed in [71 [8]] to analyze metastability through a systematic use 
of classical potential theory. Let us recall the basic notions we will need. 

For two disjoint sets A, B c Sn, the equilibrium potential, Ha,b, is the harmonic 
function, i.e. the solution of the equation 

(Lh AB )(a) = 0, a^AUB, (2.1) 

with boundary conditions 

h AB (a) = l 1 ' iiaeA . (2.2) 

' V ; [0, if a E B 

The equilibrium measure is the function 

£a,b{o) = -{Lh A)B )(a) = (Lh B>A )(a), (2.3) 

which clearly is non-vanishing only on A and B. An important formula is the 
discrete analog of the first Green's identity: Let D c S N and D c = S N \ D. Then, 
for any function /, we have 

~ £ »{°)PN{v,°')[f{v)-fW)\ 2 (2.4) 



$>(*)/(*)(£/)(*) - E vwmmw- 



METASTABILITY IN THE RFCW MODEL 



7 



In particular, for / = h AjB , we get that 

2 Yl ^{o)PN{(?,(?')[h A)B (a) - h A)B (a')f (2.5) 

= ^2fi(cr)e AtB (a) = c&p(A,B), 

(T&A 

where the right-hand side is called the capacity of the capacitor A, B. The func- 
tional appearing on the left-hand sides of these relations is called the Dirichletform 
or energy, and denoted 

* N (f) = ± ^)p N (a,a')[f(a)-f(aT- (2-6) 

As a consequence of the maximum principle, the function h Ai s is the unique mini- 
mizer of $jv with boundary conditions (12.2D . which implies the Dirichlet principle: 

cap(A,B)= inf & N (h), (2.7) 

where H A , B denotes the space of functions satisfying (I2.2D . 

Equilibrium potential and equilibrium measure have an immediate probabilistic 
interpretation, namely 

K[r A <r B] J h ^\ *°* A » B (2.8, 

An important observation is that equilibrium potentials and equilibrium measures 
also determine the Green's function. In fact (see e.g. USED), 

h A , B (a) = 2j G SN \ B {a,a')e AtB {a') (2.9) 
a'eA 

In the case then A is a single point, this relation can be solved for the Green's 
function to give 

l^{a)e a:B {a) 

This equation is perfect if the cardinality of the state space does not grow too fast. 
In our case, however, it is of limited use, since both numerator and denominator 
tend to be very close to zero for the wrong reason. However, Q2.9D remains useful. 
In particular, it gives the following representation for mean hitting times 

J2fi(cr)e AB (a)E a r B = ^ u.{a')h A>B (a'), (2.11) 

or, using definition (I1.27D 

K AB r B = — -L^ K a ') h A,B( a ')- (2-12) 

<r'£,Sjv 

From these equations we see that our main task will be to obtain precise esti- 
mates on capacities and some reasonably accurate estimates on equilibrium poten- 
tials. In previous applications B3 [71 [8l [9l ED , three main ideas were used to obtain 
such estimates: 
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(i) Upper bounds on capacities can be obtained using the Dirichlet variational 
principle with judiciously chosen test functions. 

(ii) Lower bounds were usually obtained using the monotonicity of capacities 
in the transition probabilities (Raighley's principle). In most applications, 
reduction of the network to a set of parallel 1 -dimensional chains was suf- 
ficient to get good bounds. 

(iii) The simple renewal estimate < cap(xB) was use ^ t0 bound the 
equilibrium potential through capacities again. 

These methods were sufficient in previous applications essentially because entropy 
were not an issue there. In the models at hand, entropy is important, and due to 
the absence of any symmetry, we cannot use the trick to deal with entropy by a 
mapping of the model to a low-dimensional one, as can be done in the standard 
Curie-Weiss model and in the RFCW model when the magnetic field takes only 
finitely many values Ifl5l, [H] . 

Thus we will need to improve on these ideas. In particular, we will need a 
new approach to lower bounds for capacities. This will be done by exploiting a 
dual variational representation of capacities in terms of flows, due to Berman and 
Konsowa fl2|] . Indeed, one of the main messages of this paper is to illustrate the 
power of this variational principle. 

Random path representation and lower bounds on capacities. It will be con- 
venient to think of the quantities /j l (a)p N (a, a') as conductances, c(a, a'), associated 
to the edges e = (a, a') of the graph of allowed transitions of our dynamics. This 
interpretation is justified since, due to reversibility, c(a, a') = c(a' } a) is symmetric. 

For purposes of the exposition, it will be useful to abstract from the specific 
model and to consider a general finite connected graph, (S, S) such that whenever 
e = (a, b) G S, then also — e = (b, a) G S. Let this graph be endowed with a 
symmetric function, c : S — > E + , called conductance. 

Given two disjoint subsets A,BcS define the capacity, 

cap(A,S) = - min V c(a, b) (h(b) - h(a)f . (2.13) 

2 h\ A =Q, h\ B =l 

e=(a,b)&c 

Definition 2.1. Given two disjoint sets, A, B c S, a non-negative, cycle free unit 
flow, f, from A to B is a function / : £ — > M + U{0}, such that the following 
conditions are verified: 

(i) if/(e)>0,then/(-e) = 0; 

(ii) / satisfies Kirchoff's law, i.e. for any vertex a G S \ (A U B), 

^/(6,a) = ^/(a,d); (2.14) 
b d 

(iii) 

E E /(«> 6 ) = 1 = E E /(«> 6 ); ( 2 - 15 ) 

aeA b a b£B 

(iv) any path, 7, from A to B such that /(e) > for all e G 7, is self-avoiding. 
We will denote the space of non-negative, cycle free unit flows from A to B by 
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An important example of a unit flow can be constructed from the equilibrium 
potential, h*, i.e. the unique minimizer of (I2.13D . Since h* satisfies, for any a G 

S\(AUB), 

J2 c ( a ^)(h*(b)-h*(a)) = 0, (2.16) 

6 

one verifies easily that the function, /*, defined by 



1 

cap (A, B) 



fM) = zzzn r,A a M ( h *( a ) - h*(b)) + , (2.17) 



is a non-negative unit flow from A to B. We will call /* the harmonic flow . 

The key observation is that any f G Ua,b gives rise to a lower bound on the 
capacity cap(^4, B), and that this bound becomes sharp for the harmonic flow. To 
see this we construct from / a stopped Markov chain X = (X , . . . , X r ) as follows: 
For each ae S\B define F(a) = ^ b f(a, b) . 

We define the initial distribution of our chain as (a) = F(a), for a G A, and 
zero otherwise. The transition probabilities are given by 

9 'M) = #|, (2.18) 
F{a) 

for a £ B, and the chain is stopped on arrival in B. Notice that by our choice of 
the initial distribution and in view of (I2.18D X will never visit sites a G S \ B with 

F(a) = 0. 

Thus, given a trajectory X = (a , ai, . . . , a r ) with a G A, a r G B and a £ G 
5\ {AUB) for£ = 0,...,r- 1, 

ru=o F ( a ^) 

where = (a^, a^ + i) and we use the convention 0/0 = 0. Note that, with the above 
definitions, the probability that X passes through an edge e is 

P / (e G X) = £V(AT)l {e6 * } = /(e). (2.20) 
Consequently we have a partition of unity, 



1 {/(e)>0} - 2^ J(Z\ ' ( - 2 ' 21J 



A 



We are ready now to derive our /-induced lower bound: For every function h with 

h\ A = and h\ B = 1, 



e:/(e)>0 
A eeA ''W 



METASTABILITY IN THE RFCW MODEL 



10 



As a result, interchanging the minimum and the sum, 

r-l 



c(ag, d£ + i) N L/„ \\2 

r A"=(ao,...,a r ) 

/(e) 



v , s h(ao)=0, h{ar)=l~~ J {dt, CLt+l) 



E p/ w 



X 



E 



(2.22) 



Since for the equilibrium flow, /*, 

c(e) cap(A, B) 



\-IM = I (223) 

^ r(p\ PtL-nfA RV ' ^ 



with F- f * -probability one, the bound (12.221) is sharp. 
Thus we have proven the following result from [121] : 

Proposition 2.2. Let A,BcS. Then, with the notation introduced above, 



cap(A, B) = sup E f 

fev A , B 



E 



e,X C ^ 



(2.24) 



3. Coarse graining and the mesoscopic approximation 



The problem of entropy forces us to investigate the model on a coarse grained 
scale. When the random fields take only finitely many values, this can be done by 
an exact mapping to a low-dimensional chain. Here this is not the case, but we 
can will construct a sequence of approximate mappings that in the limit allow to 
extract the exact result. 

3.1. Coarse graining. Let I denote the support of the distribution of the random 
fields. Let 1%, with t e {1, . . . , n}, be a partition of I such that, for some C < oo 
and for all £, \I e \ < C/n = e. 

Each realization of the random field {/i^a;]}^ induces a random partition of the 
set A = {1, ... , N} into subsets 

A k [cu} = {ieA: tn[u] G I k }. (3.1) 

We may introduce n order parameters 

m k [uj](a) = — ^2 a i- ( 3 - 2 ) 

ieA fc [w] 

We denote by m [u] the n-dimensional vector (mi[w], . . . , m n [uj}). In the sequel we 
will use the convention that bold symbols denote n-dimensional vectors and their 
components, while the sum of the components is denoted by the corresponding 
plain symbol, e.g. m = Y^e=i rn e- rn takes values in the set 

Y n N [u] = x% =1 {-p Ntk [u}, -p N , k [uj} + jf, . . . , p N , k [uj} - jj, p Nik [uj}} , (3.3) 

where 

P k =p N A") = ^- (3.4) 

We will denote by e^, £ = 1, . . . , n, the lattice vectors of the set T 1 ^, i.e. the vectors 
of length 2/N parallel to unit vectors. 
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Note that the random variables p N>k concentrate exponentially (in N) around 
their mean values ~K h p Nik = F h [hi e I k ] = p k . 

Notational warning: To simplify statements in the remainder of the paper, we will 
henceforth assume that all statements involving random variables on (O, J 7 , F h ) 
hold true with F h -probability one, for all but finitely many values of N. 

We may write the Hamiltonian in the form 

n 

H N [u](a) = -NE{m[u){a))+Y,H a ^^ (3 - 5) 
where E : IR n — ► R is the function 

E{x) = ]-[^2x k \ +^2h k x k , (3.6) 



=1 ieA 




k=l 



with 



hi = tt-t hi, and hi = hi - h e . (3.7) 

Note that if h t = hi for all % e A e , which is the case when h takes only finitely many 
values and the partition li is chosen suitably, then the Glauber dynamics under the 
family of functions mi is again Markovian. This fact was exploited in lfl"5~l [6D . Here 
we will consider the case where this is not the case. However, the idea behind our 
approach is to exploit that by choosing n large we can get to a situation that is 
rather close to that one. 

Let us define the equilibrium distribution of the variables m[a] 

Qf3 >N [oj](x) = pp tN [u](m[uj](a) = x) (3.8) 
- 1 r WE(x) ¥ ■* iP £?=i£ i6 a^(^-M 

where Z N [u] is the normalizing partition function. Note that with some abuse of 
notation, we will use the same symbols Qp t N, Fp^ as in Section 1 for functions 
defined on the n-dimensional variables x. Since we distinguish the vectors from 
the scalars by use of bold type, there should be no confusion possible. Similarly, 
for a mesoscopic subset A C r^[u;], we define its microscopic counterpart, 

A = S N [A] = {a eS N : m(a) G A} . (3.9) 

3.2. The landscape near critical points. We now turn to the precise computation 
of the behavior of the measures Qp^ N [uj](x) in the neighborhood of the critical 
points of Fp^ N [uS\(x). We will see that this goes very much along the lines of the 
analysis in the one-dimensional case in Section 1. 
Let us begin by writing 




Z PtN [u]Qp tN [aj](x) =exp j N(3 j - ( ] +^x t hz ] ] J[ Zp tN M(x t /pt), 

(3.10) 
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where 



Z^ N [u)(y) =E aA( exp 1(3 £ ~h m ) l^-i a ^ = E ^ V 1 ^"^} 



(3.11) 

For y G (—1,1), these can be expressed, using sharp large deviation estimates 
ILIQII, as 

4 = eM-\A e \i N ,M(y)) + (3J2) 

/ f|A,|//^M(y) 



where o(l) goes to zero as |A^| f oo. Note that as in the one-dimensional case, we 
identify functions on with their natural extensions to IR n . This means that we 
can express the right-hand side in (13.10D as 



( I N.en( x e/Pe)/Pe 

where 



ZwMQeAuK*) = II V 6X p (-^H(*)) (i + 0(1)) , 

(3.13) 



Fp )N [<J\(x) = -- I ^2x t J -^x^/ + ^^/£>//jv,/H(aj < /p/). (3.14) 

Here Jjv,i[a;](y) is the Legendre-Fenchel transform of the log-moment generating 
function, 

U N ,M{t) = i^rlnES^expU^rTi) (3.15) 

= — — - In cosh ft + 0hj \ . 

We again analyze our functions near critical points, z* , of Fp^. Equations (I3.10D - 
(I3.15D imply: if z* is a critical point, then, for ||v|| < N~ l l 2+& , 

= (1 + o(D) , (3-16) 

with 

(A(**))« = /'^ j = -1 + SwP^pjHtrtte/pi) = -1 + (3.17) 
Now, if z* is a critical point of F^ N , 

or, with z* = £J =1 z|, 

/3 (z* + ^) = r Ntt (z* t /p e ) = t\. (3.19) 
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By standard properties of Legendre-Fenchel transforms, we have that I' Ne {x) = 
U'x\{x), so that 

4/ Pi = KM** + h i)) ^TTlY] tanh(/3(z* + hi))). (3.20) 

Summing over t, we see that z* must satisfy the equation 

** = ij]tanh(/3(z* + /0), (3.21) 
ieA 

which nicely does not depend on our choice of the coarse graining (and hence on 

ri). 

Finally, using that at a critical point I'^ e (z* e /pi) = v „ 1 ^ t *y we get the explicit 
expression for the random numbers on the right hand side of A3.17D 

Xi = P Pe U^(f3(z* + h t )) = ^E ieA ,(l-tanh 2 (^ + ^)))- (3 " 22) 
The determinant of the matrix A(z*) has a simple expression of the form 

/ n \ n 

det(A(z*)) = 11^ (3 - 23) 

\ e=i A V i=i 



ieA / £=l 



e 



(1 - (3E h (1 - tanh 2 ^* + h)))) ]JX £ (1 + o(l)) 



£=1 

where o(l) J, 0, a.s., as | oo. Combing these observations, we arrive at the 
following proposition. 

Proposition 3.1. Let z* be a critical point of Qp t N- Then z* is given by (\3.20\) where 
z* is a solution of (\3.21\) . Moreover, 



, *s \f\ det(A(z*),, 

Zp,nQp,n(z ) = Vl 1 1 j/ ' (3.24) 



\/3E h (1 - tanh 2 (/3(z* + h))) - l| 
x exp ((IN + -L gin cosh {(3{z* + h,)) JJ (1 + o(l)) . 

Proq/: We only need to examine (I3.13D at a critical point 2*. The equation for the 
prefactor follows by combining (I3.12D with (I3.23D . As for the exponential term, 
Fp } N, notice that by convex duality 

I N ,i{z* e / Pi ) = t* t z* e /p t - U N)i {f e ) = f3{z* + ht)z* t /p t - U N>e (P(z* + h t )) ■ (3.25) 
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Hence (13.141) equals 

n 



i=i 

1 . „2 



2 



z* e h £ - z*z* e - hz* e + — ^ ^ In cosh + ft,*)) 



-(z*) 2 - 



^^lncosh (/?(>* + /*;)). (3.26) 



□ 



Remark. The form given in Proposition 13.11 is highly suitable for our purposes as 
the dependence on n appears only in the denominator of the prefactor. We will see 
that this is just what we need to get a formula for capacities that is independent of 
the choice of the partition of / and has a limit as n | oo. 

Eigenvalues of the Hessian. We now describe the eigenvalues of the Hessian 
matrix A(z*). 

Lemma 3.2. Let z* be a solution of the equation (\3.21 I ). Assume in addition that all 
numbers X k are distinct. Then 7 is an eigenvalue o/A(z*) if and only if it is a solution 
of the equation 

n 

<=1 # Ei £A ^ (l-tarfi a (/9(**+Ai))) 7 

Moreover, A3. 27V has at most one negative solution, and it has such a negative solution 
if and only if 

8 - 

£ (! - tanh2 (/5 (** + >**))) > 1- (3-28) 
1=1 

Remark. To analyze the case when some coincide is also not difficult. See flU] . 

Proof. To find the eigenvalues of A, just replace by — 7 in the first line of 
(IQ5D . This gives 

/ n \ n 

det (A(z*) -7))= 1-E * II& " 7), (3-29) 

V £=1 A£-7y £=1 

provided none of the A^ — 7 = 0. (13.271) is then just the demand that the first 
factor on the right of (13.291) vanishes. It is easy to see that, under the hypothesis of 
the lemma, this equation has n solutions, and that exactly one of them is negative 
under the hypothesis (I3.28D . □ 

Topology of the landscape. From the analysis of the critical points of Fp >N it 
follows that the landscape of this function is closely slaved to the one-dimensional 
landscape described in Section 1. We collect the following features: 

(i) Let m\ < z{ < m* 2 < zl < ■ ■ ■ < z* k < m* k+1 be the sequence of minima 
resp. maxima of the one-dimensional function Fp^ defined in (11.10D . Then 
to each minimum, m*, corresponds a minimum, m* of Fp jN , such that 
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Figure 1. Correspondence of one and n-dimensional landscape 

J2e=i m *i i = m i> an d two each maximum, z*, corresponds a saddle point, 
z* of F/3 lN , such that Y%=i z* A = z*. 
(ii) For any value m of the total magnetization, the function Fp tN (x) takes 
its relative minimum on the set {y : Y^Vi = m } at tne point x e M™ 
determined (coordinate-wise) by the equation 

Xi(m) = ^2 tanh (/3 (m + a + /ij)) , (3.30) 



where a = a(m) is recovered from 

1 

N 



hi = ^- ^ tanh (/5 (m + a + ^)) . (3.31) 



ieA 
Moreover, 

Fp, N (m) < Fp >N {x) < F 0tN (m) + 0(n\nN/N). (3.32) 

Remark. Note that the minimal energy curves x(-) defined by (I3.30D pass through 
the minima and saddle points, but are in general not the integral curves of the 
gradient flow connecting them. Note also that since we assume that random fields 
{hi(uj)} have bounded support, for every 5 > there exist two universal constants 

< ci < c 2 < oo, such that 

dxAm) 

ciPe < — ] < c 2 pi, (3.33) 

dm 

uniformly in N, m e [—1 + 5,1 — 5] and in I — 1, . . . , n. 



4. Upper bounds on capacities 

This and the next section are devoted to proving Theorem 11.31 In this section 
we derive upper bounds on capacities between two local minima. The procedure 
to obtain these bounds has two steps. First, we show that using test functions that 
only depend on the block variables m(a), we can always get upper bounds in terms 
of a finite dimensional Dirichlet form. Second, we produce a good test function for 
this Dirichlet form. 
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4.1. First blocking. Let us consider two sets, A, B c <Sjv, that are defined in terms 
of block variables m. This means that for some A,B C r^, A = «Sjv[A] and 
B = S N [B] . Later we will be interested in pre-images of two minima of the function 
Fp,N' We S et the obvious upper bound 

cap(A,B) = inf ^ V /i^M^p^o-') [/i(o-) - /i(o-')] 2 

heH A ,B I ~*„ 

u,u Go at 

- W f o zC ^/J,^M(<r)p(<7,<7 / ) Km(cr)) -u(m(cr'))] 2 

(T,0"'60JV 

= inf V [u(x)- u(x')} 2 V fi p>N [u]{a) V p(cr, </) 
ueGA,B A — ' z — ' ^ — ' 

a;,a;'er™ creSjvM cr'e5]v[a;'] 



inf V Qp iN [u]{x)r N {x, x') [u(x) - u(x') 



,2 



;a,b). (4.D 



with 



Here 



r N (x,x') = - — l —- — W- 2 ) 



H a ,b = {h: S n -^ [0,1] : Va G A, /i(cr) = 1 , W G £?, h(a) = 0} (4.3) 

and 

Qa.b = {u:T n N -^ [0, 1] : Vaj e A, u{x) = 1, Vaj e u(aj) = 0}. (4.4) 

4.2. Sharp upper bounds for saddle point crossings. Let now z* be a saddle 
point, i.e. a critical point of Q^n such that the matrix A(z*) has exactly one 
negative eigenvalue and that all its other eigenvalues are strictly positive. Let 
A, B be two disjoint neighborhoods of minima of Fp^ N that are connected through 
z*, i.e. A and B are strictly contained in two different connected components of 
the level set {x : Fp iN (x) < F/3 :N (z*)}, and there exists a path 7 from AtoB such 
that max £ce7 F^ N {x) = Fp fN (z*). 

To estimate such capacities it suffices to compute the capacity of some small set 
near the saddle point (see e.g. fl3]] or Q8]] for an explanation). For a given (small) 
constant p = p(N) <C 1, we define 

D N {p) = {xeT n N : \z\ -x t \<p^\<i< n}, (4.5) 

In this section we will later choose p = Cy/lnN/N, with C < 00. D N (p) is the 
hypercube in T 7 ^ centered in z* with sidelenght 2 p. For a fixed vector, v G T 7 ^, 
consider three disjoint subsets, 

W = {xeT n N :\(v,(x-z*))\<p} 
W x = {xeT n N :(v,(x-z*))<-p} 

W 2 = {xeT n N :(v,(x-z*))>p}. (4.6) 

We will compute the capacity of the Dirichlet form restricted to the set D N (p) 
with boundary conditions zero and one, respectively, on the sets W\ D D N (p) and 
W 2 n D N (p). This will be done by exhibiting an approximately harmonic function 
with these boundary conditions. Before doing this, it will however be useful to 
slightly simplify the Dirichlet form we have to work with. 
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Cleaning of the Dirichlet form. One problem we are faced with in our setting 
is that the transition rates r N (x, x') are given in a somewhat unpleasant form. At 
the same time it would be nicer to be able to replace the measure Qp,N by the 
approximation given in (I3.18D . That we are allowed to do this follows from the 
simple assertion below, that is an immediate consequence of the positivity of the 
terms in the Dirichlet form, and of the Dirichlet principle. 

Lemma 4.1. Let & N , <!> N be two Dirichlet forms defined on the same space, T, corre- 
sponding to the measure Q and transition rates r, respectively Q and f. Assume that, 
for all x, x' G r, 

l Q{x) 



Then for any sets A, B 



Q(x) 



<S, 



r(x, x 



r(x, x') 
C^(A,B 



< 5. 



<(!-*)- ■ 



Cap N (A,B) 

Proof. Note that Cap^(A, B) = ini ue g A B ® N (u) = & N (u*), and 
Cap^(A, B) = mf u& g A B $> N (u) = $ N (u*). But clearly 



(4.7) 



(4.8) 



9 22 Q{X) ^ {X > X) 7(X^) {U ^~ U ^ 



(4.9) 



x,x'er 



> \ E e(*)(i-W*,aO(i -*)(«*(*) -ttV)) 



x,x'er 



> (l-<5) 2 inf - V Q(x) f(x, x') (u(x) - u(x')) 



uGGa,b 2 

x,x'er 

= (l-6) 2 cZp n N (A,B). 
By the same token, 

®n(u*) > (l-5) 2 Cap n N (A,B). 
The claimed relation follows. 



(4.10) 

□ 



To make use of this observation, we need to control the rates r N (x, x') and the 
measure Q^n(x) in terms of suitable modified rates and measures. In fact, we see 
easily that 



Qp,n{x) = Qp, N (z*)exp 



(3N 



((x- z*),A(z*)(x- z*)) 



(4.11) 



so that, for all x G D N (p) and for some K < oo, it holds 



QffAx) 



< KNp . 



(4.12) 



For that concerns the rates, let us first define, for a G S 



A^)E{ieA fc : <r(i) = ±l} 



(4.13) 
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For all x G V%, we then have 

r N (x,x + e e ) = Qfs^ix)- 1 E u.p >N [aj](a) E K ", O (4.14) 

Notice that for all cr g S N (x), \Aj(a) \ is a constant just depending on x. Using that 
hi = h e + h i} with hi G [— e, e], we get the bounds 

r N {x,x + e t ) = ^e- 2 ^ )+h l(l + 0(e)). (4.15) 
It follows easily that, for all x G D N (p), 

r N (x,x + eA 



r N (z*, z* + eA 



1 



<(3{e + np) (4.16) 



With this in mind, we let L N be the generator of the dynamics on D N (p) with 
rates r(x,x + e £ ) = rpf(z*,z* + e^) = and r(a; + e e ,x) = ^ g^'^ an d 

thus with reversible measure Qp,n(x). For w G we write the corresponding 

Dirichlet form as 

n 

i Dn {u) = q p , n {z*) e E^ e " /37V((a, " 2 * ),A(z * )(3: ^ )) ( M ( :r )- n ( :E+ ^)) 2 - 

x&D N {p) £=1 

(4.17) 

4.3. Approximately harmonic functions for &d n - We will now describe a func- 
tion that we will show to be almost harmonic with respect to the Dirichlet form 
$d n - Define the matrix M(z*) = B with elements 

B £ , fc = ^k{z*) l)k ^T k . (4.18) 

Let v^>, i = 1, . . . ,n be the normalized eigenvectors of B, and 7, be the corre- 
sponding eigenvalues. We denote by 71 the unique negative eigenvalue of B, and 
characterize it in the following lemma. 

Lemma 4.2. Let z* he a solution of the equation (\3.21\ ) and assume in addition that 

8 N 

^ e ( x - tanh2 (p ( z * + h M) > L w- 19 ) 

i=i 

Then, z* defined through (\3.20\ ) is a saddle point and the unique negative eigenvalue 
o/B(z*) 15 the unique negative solution, tl = 71 (iV, n), of the equation 

^ ^Ei 6 A,(l-tanh(^(^ + / li )))exp(-2/5[^ + /i,] ) 

> Pt— T~- r — ^1 — = 1- (4.20) 

Jiq\ E iSA (l-tanh(/3( 2 *+^))) exp (-2p[z*+h t ] + ) 

p^EieA^l-taiih 2 ^* +/><))) 7 

Moreover, we have that 

lim lim 71 (iV, n) = 71, (4.21) 

nfoo iVToo 
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where is the unique negative solution of the equation 

1 - tanh(/3(z* + h))) exp (-2(3 [z* + h]_ 



cxp(-2/3[z*+/i] + ) _^ 



1. (4.22) 



/3(l+tanh(/3(>*+/i))) 

Proof. The particular form of the matrix B allows to obtain a simple characteriza- 
tion of all eigenvalues and eigenvectors. The eigenvalue equations can be written 
as 

n 

- \/reTkUe. + (r k X k - j)u k = 0, VI < k < n. (4.23) 
t=i 

Assume for simplicity that all r k X k take distinct values. Then there is no non-trivial 
solution of these equation with 7 = r k X k , and thus Y^i=x \[ri u t 7^ 0. Thus, 

u k = v^X^v^ (4.24) 
r k \ k - 7 

Multiplying by ■ s /r k ' and summing over k, u k is a solution if and only if 7 satisfies 
the equation 

n 

V /* = 1. (4.25) 

fc=1 r k X k - 7 



Using (14.151) and noticing that ^1 = I(p fc _ we get 

r fc = \(p k - zl) exp (-2/5 [m(a) + h k ] + } (1 + 0(e)). (4.26) 



Inserting the expressions for z* k /p k and Afc given by (13.20D and (I3.22D into 04.261) 
and substituting the result into (14.251) , we recover 04.201) . 

Since the left-hand side of G4.25D is monotone decreasing in 7 as long as 7 > 0, 
it follows that there can be at most one negative solution of this equation, and such 
a solution exists if and only if left-hand side is larger than 1 for 7 = 0. The claimed 
convergence property A4.21D follows easily. □ 

We continue our construction defining the vectors i?W by 

vf=vf/^F h (4.27) 

and the vectors by 

vf = vf^r t = n vf. (4.28) 

We will single out the vectors v = and v = v^K The important facts about 
these vectors is that 

A*W = 7iV » ) (4.29) 

and that 

(i>(0 jt> (j)) = ( j. i . (4.30) 
This implies the following non-orthogonal decomposition of the quadratic form A, 

n 

(y, Ax) =^7i(y, v®)(x,v®). (4.31) 

i=i 

A consequence of the computation in the proof of Lemma Fk2l on which we shall 
rely in the sequel, is the following: 
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Lemma 4.3. There exists a positive constant 5 > such that independently ofn, 



S < minvu < maxt)j. < — . 
k k o 



Proof. Due to our explicit computations, 

1 



2 V Pfc 



0- 



A. 



^ (1 - tanh 2 {(3{z* + fci))) 



ieA k 



-2/3\z*+h k \ 



(4.32) 



(4.33) 



Consequently, the quantities <p k = r k X k — 71 (iV, n) are bounded away from zero and 
infinity, uniformly in N, n and k = 1, . . . , n. Since by A4.27D and (I4.24D the entries 
of v are given by 



-1/2 



the assertion of the lemma follows. 
Finally, define the function / : R 



ra 

J —00 



R + by 

-W\jl\u 2 /2 du 



(4.34) 

□ 

(4.35) 



gjVjTij 
2vr 



We claim that the function 

g(x) = f((v,x)) (4.36) 
is the desired approximately harmonic function. 

Notice first, that g(x) = o(l) for all x e W x H D N (p), while g(x) = 1 - o(l) for 
all a; e W 2 H D N (p). Moreover, the following holds: 

Lemma 4.4. Let g be defined in (\4. 36 P . Then, /or all x e D N (p), there exists a 
constant c < 00 such that 



Ljv0 (oj 



< 



27riV' 



-/3JV|7i|(sd 



r e v e cp 



(4.37) 



Remark. The point of the estimate (14.37D is that it is by a factor p 2 smaller than 
what we would get for an arbitrary choice of the parameters v and 71. We will 
actually use this estimate in the proof of the lower bound. 

Proof. To simplify the notation we will assume throughout the proof that coordi- 
nates are chosen such that z* = 0. We also set A = A(z*). Using the detailed 
balance condition, we get 



r(x,x- ep) 



Q/3,n(x - ep) 
Qi3,n(x) 



rix — ep, x) 



Qp,n(x - e j) 



-n- 



(4.38) 



Moreover, from the definition of <2/3,tv and using that we are near a critical point, 
we have that 



Qp,n(x - ep) 
Qp,n(x) 



exp(-^f[(x,Ax) - {{x-e e ),A{x-e £ ))]) (4.39) 
exp (-/3 (e e , Ax) ) (l + O (N' 1 ) ) . 
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From (I4.38D and (14.39D , the generator can be written as 

n 

(Z N g)(x) = J2 r z(9(x + ee) - g(x)) (4.40) 

(l - exp { -p {et , A,)) ^ (1 + ° (Jr '») ■ 



*=1 

X 



Now we use the explicit form of g to obtain 

^(x + c^-^aj) = f((x,v)+v e /N)-f{(x,v) (4.41) 
= /'((a;, v))v t /N + v 2 N- 2 f"(x, v)/2 + v 3 e N- 3 f"'((x, v))/6 

= ^ySe-^K^/J. (1 _ ^ N(a , „)/2 + O (p 2 )) . 

In particular, we get from here that 

g{x)-g{x-e f ) = ^ (-^JV^al [(a: - c <f «) 2 - (as, «) 2 ] /2) (4.42) 
g[x + e^j 9{ x ) 

l-vM 1 \[(x,v)-v e /N]/2 + Q(p i ) 

X l-^| 7l |(a^)/2 + C>(p 2 ) 

/ or* I , u A , ^ 2 /3|7i|/2iV + 0(p 2 ) \ 
= exp^KI^,)) {^ 1 _ vemil{x ^ +0(p , ) ) 

= exp(-(3\^\v e (x,v))(l+0(p 2 )) 
Let us now insert these equations into (14.40D : 



(1 -exp{-(3(e e ,Ax) - (3\^\v e (x,v)} (l + 0(p 2 ))) . (4.43) 



mi 

£=1 

x / " 
Now 

l-exp(-/3(e*Aa:) - /3| 7 i|^K v)) (l + 0(p 2 )) 

= /?(e £ ,Aa ; )+/3|7 1 |^(a ; ,- l ;) + 0(p 2 ). (4.44) 
Using this fact, and collecting the leading order terms, we get 



(l Ng )(x) = ^m e -m^f/2 

n 

xj^nvt [((3(e i ,Ax)+(3\j 1 \vt(x,v)) +0(p 2 )] . (4.45) 
Thus we will have proved the lemma provided that 

rtvi ( (e*, Ax) -7i«<(aj, v)) = 0. (4.46) 

£=1 

But note that from (I4.31D we get that 

n 

(e<,Aa:) -7 lV< (aj,«) = ^Tjwf (aj, (4.47) 

i=2 
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Hence using that by (I4.28D r t vi = v e and that by (14.30D v is orthogonal to with 
j > 2, (14.46D follows and the lemma is proven. □ 

Having established that g is a good approximation of the equilibrium potential 
in a neighborhood of z* , we can now use it to compute a good upper bound for 
the capacity. Fix now p = C-^/lnN/N. 

Proposition 4.5. With the notation introduced above and for every n 6 N, we get 

(4.48) 

Proq/: The upper bound on cap(A,B) is inherited from the upper bound on the 
mesoscopic capacity Cap^(A, B). As for the latter, we first estimate the energy of 
the mesoscopic neighborhood D N = D N (p) of the saddle point z*. By Lemma [47TL 
this can be controlled in terms of the modified Dirichlet form <f> Dn in A4.17D . Thus, 
let g the function defined in (14.36D and choose coordinates such that z* = 0. Then 

n 

$ DN (g) = Q AJV (0) ]T Y, e ~ miX ' AX))/2r ^9^ + ee)-g(x)) 2 (4.49) 

xeD N i=i 

x£D N 1=1 

x (1 -t»//3|7i|(aj,«) + (^MnA^)) 2 



Here we used that J2e r e v j — 2~2e v 2 e = I. It remains to compute the sum over x. By 
a standard approximation of the sum by an integral we get 

^2 e -^ 7V l^l( £C '' u ) 2 e" /37V({a '' Aa,))/2 (4.50) 



1 

n I n 



i=i ) 

n In \ r 



2/3 
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Inserting (14.50D into A4.49D we see that the left-hand side of (14.49D is equal to the 
right-hand side of (I4.48D up to error terms. 

It remains to show that the contributions from the sum outside D N in the Dirich- 
let form do not contribute significantly to the capacity. To do this, we define a 
global test function g given by 

f 0, x G W\ 
g(x) = I 1, xeW 2 (4.51) 
[ g(x), x G Wq 

Clearly, the only non-zero contributions to the Dirichlet form &n{9) come from 
W = W U dWo, where dW denotes the boundary of W . Let us thus consider the 
sets W™ = W n D N and W° ut = W nD c N (see Figure 03]). We denote by ^ wm {g) 
the Dirichlet form of g restricted to Wq 1 and to the part of its boundary contained 
in D N , i.e. to W\ n n D N , and by <&^ out (g) the Dirichlet form of g restricted to W° Q ut . 
With this notation, we have 



$ N (g) = ® l l t (g) + $^r S ut(g) (4.52) 
®l t (g) - ($l t (g) - $l t (g))) (l + O (V^n/n)) + ^ r (g). 



The first term in (I4.52D satisfies trivially the bound 

$ D ' N (g) < ® l U(g) < *dM, (4-53) 



where D' N = D N (p') is defined as in (l4~53l) but with constant p' = C'y/lnN/N 
such that D' N c W™. Performing the same computations as in (I4.49D and (14.501) 
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it is easy to show that §D' N {g) = ®D N {g){^ + o(l)), and then from (I4.52D it follows 
that 

*i*>(g) = *DM(l-o(l)). (4.54) 
Consider now the second term in A4.52D . Since g = g on W , we get 

n 

®w i M-®wUg)= E E QipYi [(9(x + e t ) - g{x)f - g(x) 2 ] 

xedw^nWi e=i 

n 

+ E E^^ + ^-rf-^-rf]- C4.55) 

where we also used that the function g has boundary conditions zero and one 
respectively on W\ and W 2 . By symmetry let us just consider the first sum in the 
r.h.s. of (14351) . For x e dW l Q n n W x it holds that (x, v) < -p = -Cy/\nN/N, and 
hence 

g ( x f < 1 — e -m^\p\ ( 4.5 6 ) 

~ y/2irP\*ii\Cy/hrN 

Using this bound together with inequality (14.41D to control [g{x + e ( ) — g(x)) 2 , we 
get 

E E GO")** to (* + ^ - ^)) 2 - ^*) 2 ] 



cceaw^nWi \ v / 

< Q (o)^-^ y c -/W((M-))/a ( ! + (4.57) 

scsaw * n rWi v v v 

for some constant c independent on AT. The sum over a? e dW£ n fi Wi in the last 
term can then be computed as in G4.50D . However, in this case the integration runs 
over the (n — 1) -dimensional hyperplane orthogonal to v and thus we have 



e -/3N(( X ,Ax))/2 



E 

kr ( " 

7 V=2 
n— 1 / n 



1=2 J J 

< fjY' 1 lf[Vnj e~^ 2 / 2 J d^ye-^&U^*) 2 /*) 
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Inserting (I4.58D in A4.57D , and comparing the result with $ Djv (g), we get that the 
l.h.s of A4.57D is bounded as 

i 1 + C \^n) ^ e ~ mhllp2/2 ®»M = o(N- K )$ DN (g), (4.59) 

with K = ffmjgzl ^ w hich is positive if C is large enough. A similar bound can be 
obtained for the second sum in (I4.55D , so that we finally get 

$l t (g) - $ l l t (g)\ < o(N- K )$ DN (g). (4.60) 

The last term to analyze is the Dirichlet form §^ out (g). But it is easy to realize that 

this is negligible with respect to the leading term. Indeed, since for all x 6 D C N it 
holds that Fp^ N (x) > Fp jN (z*) + K'hiN/N, for some positive K' < oo depending 
on C, we get 

*5y*(fl0 < Z^ N e-P NF ?M**) N -(K'-n) = o(iV-^")$ Djv (^). (4.61) 

From (14321) and the estimates given in (14341) , (14391) and (I4.61D , we get that 
®N(g) = ®D N (g)0- + o(l)) provides the claimed upper bound. □ 

Combining this proposition with Proposition 13.11 yields, after some computa- 
tions, the following more explicit representation of the upper bound. 

Corollary 4.6. With the same notation of Proposition \4.5\ 

Zp iN c&p(A, B) < — — , (4.5ZJ 

^ v y flNKh (1 - tanh 2 (/? (** + h))) - 1 

where 71 is defined through Eq. H4.22\) . 
Proof. First, we want to show that 



det(A(z*))|=mr, H% (4.63) 




To see this, note that 

B = RA(z*)R, 

where R is the diagonal matrix with elements Re :k = 5k,ey/re- Thus 

n n 

JJ|^| = |det(B)| = \det{RA{z*)R) \ = | det(A(z*))| det(i? 2 ) = |det(A(z*))| \\r t 

e=i i=i 

(4.64) 

as desired. Substituting in (14.481) the expression of Qp,n(z*) given in Proposition 
(13. ID . and after the cancellation due to (14.631) . we obtain an upper bound which is 
almost in the form we want. The only n-dependent quantity is the eigenvalue 71 of 
the matrix B. Taking the limit of n — > 00 and using the second part of Lemma [4721 
we recover the assertion (I4.62D of the corollary. □ 

This corollary concludes the first part of the proof of Theorem 11.31 The second 
part, namely the construction of a matching lower bound, will be discussed in the 
next section. 
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5. Lower bounds on capacities 

In this section we will exploit the variational principle form Proposition 12.241 to 
derive lower bounds on capacities. Our task is to construct a suitable non-negative 
unit flow. This will be done in two steps. First we construct a good flow for the 
coarse grained Dirichlet form in the mesoscopic variables and then we use this to 
construct a flow on the microscopic variables. 

5.1. Mesoscopic lower bound: The strategy. Let A and B be mesoscopic neigh- 
borhoods of two minima tua and m B of Fp jN , exactly as in the preceding section, 
and let z* be the highest critical point of Fp jN which lies between m A and m B . 
It would be convenient to pretend that m A , z*, m B e T^: In general we should 
substitute critical points by their closest approximations on the latter grid, but the 
proofs will not be sensitive to the corresponding corrections. Recall that the energy 
landscape around z* has been described in Subsection 3.2. 

Recall that the mesoscopic capacity, Cap^(A, B), is defined in (14.11) . We will 
construct a unit flow, f a,b, from A to B of the form 

fA,B[X,X)= — -pzr 4>a,b(x,X), (5.1) 

such that the associated Markov chain, MP^' B , Xa,b), satisfies 

P^' s I Yl ^AB(e) = 1 + o(l) 1=1- o(l). (5.2) 

In view of the general lower bound (I2.22D . Eq. (15.21) implies that the mesoscopic 
capacities satisfy 

Cap W A,B)>E!r{ £ Q J: ( tl(e) } d-(l)). C5.3) 

which is the lower bound we want to achieve on the mesoscopic level. 

We shall channel all of the flow f A ,B through a certain (mesoscopic) neighbor- 
hood G N of z* . Namely, our global flow, f a ,b, in (15. ID will consist of three (match- 
ing) parts, f A , f an d fs> where f A will be a flow from A to dG N , f will be a flow 
through Gn, and f B will be a flow from OGn to B. We will recover (15. 2D as a 
consequence of the three estimates 

PUJ>(e) = l + o(l)] =l-o(l), (5.4) 



whereas, 



( Me) =o(l)j =l-o(l) and P& ( £ B (e) = o(l) ) = 1 - o(l). 

(5.5) 

The construction of f through G N will be by far the most difficult part. It will rely 
crucially on Lemma \4A\ 
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5.2. Neighborhood G N - We chose again mesoscopic coordinates in such a way 
that z* = 0. Set p = N~ l l 2+5 and fix a (small) positive number, v > 0. Define 

G N = G N (p, v) = D N (p) n {x : (x,v) G (-vp,vp)} , (5.6) 

where i) = flW is defined in (14.281) . and £>at is the same as in (I4.5D . Note that in 
view of the discussion in Section 4, within the region G N we may work with the 
modified quantities, Q^n and rf, I — 1, . . . , n, defined in (14.111) and (14.171) . 

The boundary dG N of G N consists of three disjoint pieces, dG N = OaG n U 
d B G N U d r G N , where 

9aGn = {x E dG N : (x, v) < —up} and d B G N = {x G dG N : (x, v) > vp\ . 

(5.7) 

We choose v in (I5.6D to be so small that there exists if > 0, such that 

F PjN {x) > Ff3 iN (0) + Kp 2 , (5.8) 

uniformly over the remaining part of the boundary x G d r G^. 

Let <jr be the approximately harmonic function defined in A4.36D and (I4.51D . Pro- 
ceeding along the lines of A4.49D and A4.50D we infer that, 

$^($(1 + 0(1))= Q?A X ) n(g{x + e,)-g{x)f, (5.9) 

xeG N ud A G N tel GN {x) 

where I Gn (x) = {I : x + G Gn}- For functions, <\>, on oriented edges, (x,x + 
e^), of D N , we use the notation 4>e(x) = <f>(x, x + e f ), and set 

Fi[<j)]{x) = Qf3,N(x)r(,(j) e (x), 

n 

d?[<f>]{x) = (H<f>Kx) - H<f>Kx - e*)) • 
e=i 

In particular, the left hand side of A4.37D can be written as |djF[Vg]|/ Qp )N (x). 

Let us sum by parts in (15.91) . By (15. 8D the contribution coming from 8 t Gn is 
negligible and, consequently, we have, up to a factor of order (1 + o(l)), 

J2 g(x)dr[Vg\(x) + J2 E ftWrtix)- (5-10) 

xeG N xed A G N ^Ig n (a:) 

Furthermore, comparison between the claim of Lemma 14.41 and (14.49D (recall that 
p2 _ jy25-i ^ jv -1 / 2 ) shows that the first term above is also negligible with respect 
to $iv(fiO- Hence, 

^(SO (l + o(l))= £ £ ^i[Vflf|(a!). (5.11) 

5.3. Flow through Gjy The relation (15. IIP is the starting point for our construc- 
tion of a unit flow of the form 

U(x) = -^—T^x) (5.12) 

through Gn- Above c=l + o(l)isa normalization constant. Let us fix < u <C v 
small enough and define, 



N — 





(x, v)v 




[■■■ 


X 11-112 






\\ V \\ 





(5.13) 
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Thus, G° N is a narrow tube along the principal ^-direction f Figure |5"T3T) . We want 

to construct <f) in (I5.12D such that the following properties holds: 

PI: f is confined to G N , it runs from OaGn to d B G N and it is a unit flow. That is, 

\/x G G N ,cLF[<f>](x) = and ^( x ) = L ( 5 - 14) 

xed A G N £&I Gn (x) 

P2: is a small distortion of Vg inside G° N , 

Mx) = Vi9(x) (l + o(l)), (5.15) 
uniformly in x G G° N and £ = 1, . . . , n. 

P3: The flow f is negligible outside G° N in the following sense: For some k > 0, 

max max f Ax) < ——. (5.16) 

*eG N \G% i HV ' - N- 

Once we are able to construct f which satisfies P1-P3 above, the associated Markov 

chain ^Pjy, Xj obviously satisfies (15. 4D . 

The most natural candidate for <p would seem to be Vg. However, since g is not 
strictly harmonic, this choice does not satisfies Kirchoff 's law, and we would need 
to correct this by adding a (hopefully) small perturbation, which in principle can 
be constructed recursively. It turns out, however, to be more convenient to use as 
a starting choice 



<pf\x) = v^^Aexp (-PN\^\(x,vf/2) , (5.17) 

which, by (14.411) , satisfies 

<j>f\x) = (g(x + e t ) - g{x)) (1 + O(p)) , (5.18) 

uniformly in G N . Notice that, by (I5.12D . this choice corresponds to the Markov 
chain with transition probabilities 

q{x, x + e e ) = =?V(1 + o(l)) = q t {l + o(l)). (5.19) 

From (I3.16D and the decomposition (I4.31D we see that 



1 -J^-H<t> { ° ) ) = r.vJ^e^i-^imx^r + ix^x))) 

2nN 



In particular, there exists a constant xi > such that 

^M<e*p(- Xl A-), (5.20, 

uniformly in x G G N \ G° N and I — 1, ... ,n. 

Next, by inspection of the proof of Lemma [4741 we see that there exists X2> such 
that, 

\d^}(x)\ < X2 P 2 H^ 0) \(x), (5.21) 
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uniformly in i e G N and £ = 1, . . . ,n. Notice that we are relying on the strict 
uniform (in n) positivity of the entries v £ , as stated in Lemma |4T3l 

Truncation of Vg, confinement of f and property PI. Let C + be the positive 
cone spanned by the axis directions ei,...,e n . Note that the vector v lies in the 
interior of C + . Define (see Figure 15. 3D 

G N = int (d B G% - C+) n G N and G 2 N = {d A G l N + C+) n G N . (5.22) 

We assume that the constants v and u in the definition of Gn and, respectively, in 
the definition of G° N are tuned in such a way that G 2 N fi d r GN = 0- Let be the 




Figure 3. Narrow tube G° M and sets G\ r and G\ 



restriction of 0(°) to G N , 

$\x) = $\x)l {aefnf y (5.23) 

Now we turn to the construction of the full flow. To this end we start by setting 
the values of 4>t on 8aG n equal to 0( O ) if ^ e Ig n ( x ) an d zero otherwise. By (15. 1 ID 
and the bound (15.20D . the second of the relations in G5. 14D is satisfied. 

In order to satisfy Kirchoff's law inside Gn, we write <fi as = 0^°^ + u with u 
satisfying the recursion, 

n n 

£V*[i*](aj) = X^MO* " " d^[0 (O) ]N. (5.24) 

Since 0^°^ = on G N \ G l N , we may trivially take u = on \ G^r and then solve 
A5.24D on G 2 N using the latter as an insulated boundary condition on dG 2 N n G^. 

Interpolation of the flow inside G 2 N . We first solve A5.24D inside G N . By con- 
struction, if x G G N then x — e e <E G l N U cUGjy, for every £ = 1, . . . , n. Accordingly, 
let us slice G l N into layers L fe as follows: Set 

L = d A G N , (5.25) 
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and, for k — 0,1 



* 5***3 

k 



L fc+ i = | x G G;v : ^ - e £ G [J Lj for all £ = 1, ... ,nj . (5.26) 
Since all entries of v are positive, there exists X3 = 03(71) and M < Xz/ P> such that 

M 

i=o 

Now define recursively, for each a; G L fc+i , 

^i[w](aj) = qe (^J2^[u](x - e 3 ) - dT[^](x)j , (5.28) 

where the probability distribution, q t , . . . , q n , is defined as in (I5.19D . Obviously, this 
produces a solution of A5.24D . The particular choice of the constants q t in A5.19D 
leads to a rather miraculous looking cancellation we will encounter below. 
Properties P2 and P3. We now prove recursively a bound on u that will imply that 
Properties P2 and P3 hold. Let Ck be constants such that, for all y G 

\Ft[u](y)\ <c k p 2 F e [Vg\(y)- (5-29) 
Then, for x G L fc+1 , we get by construction (I5.28D and in view of (15.21D that 

W^l < yim^-e^l 2 

^[0 (O) ](^) " , ^[0 (O) K*) 
By our choice of in (T5T23D . 



%xp J /3JV^^(a J| «W)(c i ,«W) 1 (1 + (1/N)) 



t=2 



g J - + 2/?(e J -,t;)E^ 2 (e J -,« (0 )(a ; ,t;W) (l + ^ 2 



(l + 0(p 2 ))* 



However, for each z = 2, 



J2(e„v)(e„v^)=0. (5.32) 



Therefore, with the choice q t = ^ Ve . (1 + o(l)), we get 

F0*)]( x - e ; 



V ^3*-^ = 1 + (p% (5.33) 
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uniformly in x e G X N and / = 1, . . . , n. Thus, the coefficients q. satisfy the recursive 
bound 

Ck+i < c k (1 + 0(p 2 )) + X2P 2 , (5.34) 
with c = 0. Consequently, there exists a constant, c, such that 

c k < kp 2 ce kcp \ (5.35) 

and hence, since M < Xs/p> °m = O(p). As a result, we have constructed u on 
such that 

\T t [u]{x)\=0{p)T t [Vg\{x), (5.36) 

uniformly in as e Gjy and £ = 1, . . . , n. In particular, (I5.15D holds uniformly in 
x G G l N and hence, by (15.20D . P3 is satisfied on G X N \ G° N . Moreover, since by 
construction = on Gn \ G 2 N , P3 is trivially satisfied in the latter domain. Hence 
both P2 and P3 hold on G X N U (G N \ G 2 N ). 

It remains to reconstruct u on G 2 N \ G l N . Since we have truncated Vg outside G l N , 
Kirchoff's equation 05.241) , for x e G 2 N \G l N , takes the form jF[w](a:) = 0. Therefore, 
whatever we do in order to reconstruct 0, the total flow through G 2 N \ G l N equals 

i±g> £ X>*](*)V«^}' ^ 

By (15361) and ([57201) . the latter is of the order O ^""e^ 1 ^) . Thus, P3 is estab- 
lished. 

5.4. Flows from A to d A G N and from d B G N to B. Let f be the unit flow through 
Gn constructed above. We need to construct a flow 

f A (aJ, y = 1 + o 1 p ' A (aj, y) (5.38) 

from A to O a Gn and, respectively, a flow 

f B (aj, y) = (1 + o(l)) —-^7 4>B[x,y) (5.39) 

from ObGn to £?, such that (15. 5D holds and, of course, such that the concatenation 
f a,b = {f a, f , fs} complies with Kirchoff's law. We shall work out only the f A-case, 
the f B -case is completely analogous. 

The expressions for $jv(sO and Qp yN (x) appear on the right-hand sides of Q4.48D 
and (I3.13D . For the rest we need only rough bounds: There exists a constant 
L = L(n), such that we are able to rewrite (I5.38D as, 

<j> A (x, y) = ( 1 + °( 1 ))^(9)U(x,y) ^/^-^M^-^M')), ( 5 .40) 
Qf3,N{x)r N {x,y) 

This would imply a uniform stretched exponentially small upper bound on 4> A at 
points x which are mesoscopically away from z* in the direction of VF^n, for 
example for x satisfying 

F PiN (z*) - F p>N {x) > cN 25 - 1 . (5.41) 

With the above discussion in mind let us try to construct f A in such a way that 
it charges only bonds (x, y) for which Q5.41D is satisfied. Actually we shall do 
much better and give a more or less explicit construction of the part of f A which 
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flows through G° N : Namely, with each point x G OaG° n we shall associate a nearest 
neighbor path 7^ = (-f x (-k A (x)), . . . ,7^(0)) on such that (T5~4TT) holds for all 
y G and, 

7 a '(-A; A (a;)) G A, 7 a; (0) = a; and m(y"(- + 1)) = m(y"(-)) + 2/JV. (5.42) 
The flow from A to 8aG° n will be then defined as 

U(e)= l {ee7 ^} fc(aj). (5.43) 

By construction f A above satisfies the Kirchoff's law and matches with the flow f 
through G N on 8aG° n . Strictly speaking, we should also specify how one extends 
f on the remaining part 8aG n \ 8aG° n . But this is irrelevant: Whatever we do the 
F£' B -probability of passing through 8aGn \ 8aG° n is equal to 

= o(l). (5.44) 

x ed A G N \d A G% i 

It remains, therefore, to construct the family of paths {<y x } such that Q5.41D holds. 
Each such path 7^ will be constructed as a concatenation 7^ =711^. 

STEP 1 Construction of 7. Pick 5 such that 5 — 1 < vtla = m(rriA) and consider 
the part x[8 — 1,2*] of the minimal energy curve as described in (13.30D . Let 7 
be a nearest neighbor r^-approximation of x[5 — 1, z*], which in addition satisfies 
m(j(-+l)) = m(j(-))+2/N. Since by (I3.33D the curve x[8— 1, z*] is coordinate-wise 
increasing, the Hausdorff distance between 7 and x[5 — 1, 2;*] is at most 2^/n/N. 
Let cc A be the first point where 7 hits the set D N (p), and let tx A be the last point 
where 7 hits A (we assume now that the neighborhood A is sufficiently large so 
that u A is well denned). Then 7 is just the portion of 7 from u A to tc A . 

STEP 2 Construction of 77^. At this stage we assume that the parameter v in (15.61) 
is so small that Gn lies deeply inside D N (p). In particular, we may assume that 

F/3,n{x A ) < min {F^ N (x) : x G d A G° N ] , 

and, in view of (I3.33D . we may also assume that 

xf < x e \/x G <9aG5v an d ^ = 1, • • • , n. (5.45) 

Therefore, a? — x A has strictly positive entries and, as it now follows from C4.29D , 

(Av, x - x A ) = (v,x- x A ) > 0. 

By construction G° N is a small tube in the direction of v. Accordingly, we may as- 
sume that (Ax, x — x A ) > uniformly on 8aG° n . But this means that the function 

t : [0, 1] 1 — > (A{x A + t{x - x A ), (x A + t(x - x A )) 

is strictly increasing. Therefore, is, up to negligible corrections, increasing on 
the straight line segment, [x A , x] c IR n which connects x A and x. Then, our target 
path r] x is a nearest neighbor r^-approximation of [x A , x] which runs from x A to 
x . In view of the preceeding discussion it is possible to prepare r/ x in such a way 
that Ff3 >N (z*) - Fp >N (-) > cN 26 ' 1 along rf. Moreover, by (fB~45T) it is possible to 
ensure that the total magnetization is increasing along if . 

This concludes the construction of a flow f a,b satisfying [53] □ 
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In the sequel we shall index vertices of 7^ = 7 U rf as, 

1 x = (T(-k A ),...l x m- (5-46) 

Since, 

Fp, N {v) < F ^(z*) - Cl (y-z*, vf , (5.47) 
for every y lying on the minimal energy curve x[8 — 1, z*\ and since the Hessian of 

Fp }N is uniformly bounded on x[5 — 1, z*], we conclude that if u is chosen small 
enough, then there exists c 2 > such that 

Fp, N (j X (-)) < FpAz*) ~ c 2 M") ~ **, , (5-48) 
uniformly in x G d A G%. Finally, since the entries of v are uniformly strictly posi- 
tive, it follows from (I5.48D that, 

F AiV (7 ai (-A ; )) < F^**) - c 3 ^ (5.49) 

uniformly in a; G 9a and fc G {0, . . . , /ca(^)}- 

5.5. Lower bound on cap(A,B) via microscopic flows. Recall that A and 
are mesoscopic neighborhoods of two minima of F^n> z * 1S the corresponding 
saddle point, and A = S N [A], B = S N [B] are the microscopic counterparts of A 
and B. Let f AB = {f A , f, f B } be the mesoscopic flow from A to B constructed 
above. In this section we are going to construct a subordinate microscopic flow, 
Ja,b, from A to B. In the sequel, given a microscopic bond, b = (a, a'), we use 
e(b) = (m(a),m(a')) for its mesoscopic pre-image. Our subordinate flow will 
satisfy 

fA,*(e)= /a,b(&). (5.50) 

fe:e(6)=e 

In fact, we are going to employ a much more stringent notion of subordination on 
the level of induced Markov chains: Let us label the realizations of the mesoscopic 
chain X AB as x = (x-t A , . . . , X£ B ), in such a way that x_£ A G A, xt B G B, and 
m(aj ) = m(z*). If e is a mesoscopic bond, we write e G x if e = (a^,a^ + i) for 
some £ = — £a, . . . , — 1. To each path, a?, of positive probability, we associate a 
subordinate microscopic unit flow, f-, such that 

/^(6) > if and only if e(b) e x. (5.51) 

Then the total microscopic flow, f A ,B, can be decomposed as 

f A)B = F n' B ( x Ab = x) (5.52) 

a; 

Evidently (I5.50D is satisfied: By construction, 

/-(&) = 1 for every x and each e g x. (5.53) 

b:e(b)=e 

On the other hand, f AB (e) = £^' s (^,b = a) l{ee £ }. 
Therefore, (15.52D gives rise to the following decomposition of unity, 

Mf A , B (»o} = 2. 2. f AB (e(&))/*(&) ' (5 - 54) 

xBe(b) q_3b >A,H\ \ U 
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where (P-, S) is the microscopic Markov chain from A to B which is associated to 
the flow /— . 

Consequently, our general lower bound A2.24D implies that 



cap(A, B) > 



> 



E- 



£=-e A 

u,b(^,^+i)/-( 



P<P,n(P£)Pn{P£, at+i) 
We need to recover $7v(g) from the latter expression. In view of (15.1D . write 

Ia,b( x ^ x £+i)fH a £, ve+i) <pA,B( x e, x e+i) 



-1 
-l 

(5.55) 



Q(3, N (x e )r N (x e ,x E+1 )p 



(5.56) 



°~£+l, 



Since we prove lower bounds, we may restrict attention to a subset of good real- 
izations x of the mesoscopic chain X A B whose P|^ B -probability is close to one. 
In particular, A5.4D and (I5.5D insure that the first term in the above product is pre- 
cisely what we need. The remaining effort, therefore, is to find a judicious choice 
of /- such that the second factor in (15.561) is close to one. To this end we need 
some additional notation: Given a mesoscopic trajectory x = (x_£ A , . . . , xg B ), de- 
fine k = k(£) as the direction of the increment of £-th jump. That is, x t+ i = xt + e k . 
On the microscopic level such a transition corresponds to a flip of a spin from the 
A k slot. Thus, recalling the notation A^(cr) = {i G A fc : a(i) — ±1}, we have that, 
if a e G S N [x e ] and a £+1 G S N [x i+1 ], then a i+1 = Ofa? for some i G A~ {e) (a e ). By 
our choice of transition probabilities, p N , and their mesoscopic counterparts, r N , 
in (EL2D, 

r N (x £ , xe +1 ) 



p N (a e ,a i+ i) 



A 



k(£) 



Pi 



;i + 0(e)) 



(5.57) 



uniformly in £ and in all pairs of neighbors o^ +1 . Note that the cardinality, 



A 



, is the same for all at G Sn [xt] 



k(£)V 

For x G T^r, define the canonical measure, 



l{<7e5 J vN}^/3,Ar(c r ) 



Qp,n( x ) 



The second term in (15.561) is equal to 

f~ipii <Tt+l) 



H\nPi) 



(a t ) ■ 1/ A fcW (<r/) 



;i + 0(e)). 



(5.58) 



(5.59) 



If the magnetic fields, h, were constant on each set I k , then we could chose the 

flow f-(ae,at + x) = l^ N {^t) ■ 1/ A~,^(ai) , and consequently we would be done. 

In the general case of continuous distribution of h, this is not the case. However, 
since the fluctuations of h are bounded by 1/n, we can hope to construct /- in such 
a way that the ratio in (I5.59D is kept very close to one. 
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Construction of /-. We construct now a Markov chain, P-, on microscopic trajec- 
tories, £ = {(T , . . . , cr eB }, iromS[x ] to B, such that a e G <S[ccg], for all t = 0, ... , £ B - 
The microscopic flow, is then defined through the identity P- (b G S) = f-(b). 

The construction of a microscopic flow from A to S[x ] is completely similar (it 
is just the reversal of the above) and we will omit it. 

We now construct P-. 
STEP 1. Marginal distributions: For each t = 0, . . . ,£ B we use vj to denote the 
marginal distribution of og under P— . The measures vf are concentrated on <S[a^]. 
The initial measure, vf, is just the canonical measure The measures vf +l are 
then defined through the recursive equations 

^1(^+1)= ^IVM^^+i)- (5.60) 

a e £S[x e ] 

STEP 2. Transition probabilities. The transition probabilities, qe(a e , cr<?+i), in (15.601) 
are defined in the following way: As we have already remarked, all the microscopic 
jumps are of the form a e 1— > Of^e, for some j G A k ^(a), where flips the j-th spin 
from —1 to 1. For such a flip define 

2/3/1, 

*("«.«W = ~ PA- ( S - 6 « 

Z^ieA fc fo) c 

Then the microscopic flow through an admissible bound, b = (cr^, at + \), is equal to 
f*{a e , a t +i) = P~ (fe G £) = vl{G t )q t {a h a l+l ) = (1 + 0(e)) . (5.62) 



A fc(/)(°"<) 



Consequently, the expression in Q5.59D , and hence the second term in (15.561) . is 
equal to 



(1 + 0(e)) = ^ e (a e )(l + 0(e)). (5.63) 



Main result. We claim that there exists a set, T A B , of good mesoscopic trajectories 
from A to B, such that 

P^' s (Xa,b e T AB ) = 1 - o(l), (5.64) 
and, uniformly in x G T A B , 

E ~ [ X] M^^Bixt, x e+1 ) \ < 1 + 0(e). (5.65) 
This will imply that, 

cap(A, 5) > $ N (g) (1 - 0(e)) , (5.66) 

which is the lower bound necessary to prove Theorem ll.3l 

The rest of the Section is devoted to the proof of (I5.65D . First of all we derive 
recursive estimates on for a given realization, x, of the mesoscopic chain. After 
that it will be obvious how to define Ta b- 
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5.6. Propagation of errors along microscopic paths. Let x be given. Notice that 
/jLp e N is the product measure, 



»?,n = ®»?!£\ (5.67) 

3=1 

where is the corresponding canonical measure on the mesoscopic slot = 

{— 1, l} Aj . On the other hand, according to (I5.61D . the big microscopic chain E 
splits into a direct product of n small microscopic chains, E^, . . . , E( n ), which in- 
dependently evolve on . . . ,S^'. Thus, k{£) = k means that the £-th step of 
the mesoscopic chain induces a step of the k-th small microscopic chain E^. Let 
Ti [£],... , r n [£] be the numbers of steps performed by each of the small microscopic 
chains after £ steps of the mesoscopic chain or, equivalently after i steps of the big 
microscopic chain E. Then the corrector, in (15.631) equals 

n 
3=1 

where crp is the projection of a e on <Sjy . Therefore we are left with two separate 
tasks: On the microscopic level we need to control the propagation of errors along 
small chains and, on the mesoscopic level, we need to control the statistics of 
Ti [£],■■■ , r n [£\. The latter task is related to characterizing the set, T A B , of good 
mesoscopic trajectories and it is relegated to Subsection 15.71 

Small microscopic chains. It would be convenient to study the propagation of 
errors along small microscopic chains in the following slightly more general con- 
text: Fix 1 « M G N and < e < 1. Let gx,...,g M £ [-1,1]. Consider spin 
configurations, £ e S M = {— 1, 1} , with product weights 

w (£) = e e ^^ (i) . (5.69) 

As before, let A ± (£) = {i : £(z) = ±1}. Define layers of fixed magnetization, 
S M [K] = {£ G S M ■ |A + (OI = K}. Finally, fix S , ^ e (0, 1), such that S < S x . 

Set K = [5 M\ and r = [(5 1 - S )M\. We consider a Markov chain, 5 = 
{H , Si, ... , 5 r } on S M , such that E T e S M [K +t) = S T M for r = 0, 1, . . . , r. Let [i T 
be the canonical measure, 

MO = f 1 ^- (5-70) 

T 

We take u = fi as the initial distribution of S and, following (15.611) . we define 
transition rates, 



*r(Zr,0}Zr) = ^ ^ (5-71) 



We denote by P the law of this Markov chain and let v r be the distribution of S T 
(which is concentrated on S\.^), that is, is T (0 = P(2 T = £). The propagation of 
errors along paths of our chain is then quantified in terms of tp T (-) = z/ T (-)//i T (-). 
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Proposition 5.1. For every r = 1, . . . , r and each £ 6 <S^ define 



M 



i=l 



i=l 



(5.72) 

Then there exists c = c(5 , 5±) such that the following holds: For any trajectory, £ = 
(£o, • • • , £,r), of positive probability under P, it holds that 

A 



3 cer 2 /M 



(5.73) 



for all t = 0,1, ... ,r. 

Proof By construction, xp = 1. Let £ T+ i e Sj^" 1 . Since v T satisfies the recursion 



&V+i(£t+i) = ^(0^ £ T+ i)qv(0j £t+i,6t+i' 

jeA+(? T+ i) 

it follows that ip T satisfies 



(5.74) 



Vv+i(6m 



E 

j6A+(f T+ i) 



A*T+l(6 



T+lj 



EM^j fr+l)^^ 6r+l,£r+l) . N 
77 \ WrKPj £r+l, 

/^r+l(?r+lj 



jeA+(£ T+1 ) 

By our choice of transition probabilities in (15.711) , 



^{Qj Cr+l)gr(^- £r+l, 6r+l) _ Z T+ l 
/Vfl(£r+l) 



E 



(5.75) 



ieA-(e-? T+1 ) 



Recalling that |A + (£ r )| = |A+| = K + r does not depend on the particular value of 

fr, 



j6A+(£) 



iE^o-p^f E 



i6A-(0 



We conclude that the right hand side of (I5.75D equals 



eA-(-) 



|A + (£r+i)| £ ieA -(^ +l) e2q71 |A + (^+i)l B t (0JUiY 



(5.76) 



As a result, 



Vv+i(6 



T+l; 



|A+(€ 



jeA+(£ T+1 ) ' w 

Iterating the above procedure we arrive to the following conclusion: Consider the 
set, Z>(£ T+ i), of all paths, £ = (£ , • • • , £r, £ T +i), of positive probability from S° M to 
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S^f 1 to £ r+1 . The number, D T+ i = |P(£ T+ i)|, of such paths does not depend on 
£ T +i. Then, since ip = 1, 



We claim that 



VV+l(&H 



An 



(5.78) 



1 + 



O(e) 



A s -1 



M ) B 8 _i(e s -i 



(5.79) 



uniformly in all the quantities under consideration. Once A5.79D is verified, 



Vv(£r) < e 



0(e)r 2 /Af 



max 



L«o(€o)J 



(5.80) 



where for £ e 5^, the relation £ ~ £ r means that there is a path of positive 
probability from £ to £ T . But all such £ 's differ at most in 2r coordinates. It is then 
straightforward to see that if £ ~ £ r and £q ~ £ T , then 



< e 



0(e)r/M 



(5.81) 



and (157731) follows. 

It remains to prove (15.79D . Let £ e 5| f and <f = e S^T 1 . Notice, first of all, 
that 

B a _i(0 - = e 26 * = 1 + 0(e). (5.82) 

Similarly, 

M 

i=i 

M 
i=l 

By usual local limit results for independent Bernoulli variables, 

^_x(i G A~) - fi a (i EA-) = (J-) , (5.83) 

uniformly in s — 1, . . . , r — 1 and z = 1, . . . , M. Hence, A s -x — A s = 1 + 0(e). 
Finally, both A a -i and B s _!(£') are (uniformly ) 0{M), whereas, 

M 



A s -i - B s -i(0 = J2 - X ) e A ~) - i^eA-K')}} = °( e ) M - ( 5 - 84 ) 

Hence, 



i=l 



A-i - 1 + 0(e) 



A-i 



B s (0 S s _ 1 (e')-l + 0(e) B s _!(e') 



1 + 



0(e) 
M 



which is (157791) . 



(5.85) 

□ 
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Back to the big microscopic chain. Going back to 05.681) we infer that the correc- 
tor of the big chain S satisfies the following upper bound: Let a — (cr , cr a , . . . ) be a 
trajectory of £ (as sampled from P^). Then, for every £ — 0, 1, . . . , £ B — 1, 



Mi 



n 



.4 



(i) 
o 



1*3 W 



where Mj = \Aj\ = pjN, 

A^ = J2^^f(ieAj), and = £ 



(5.86) 



(5.87) 



■ieAj 



Of course, A^ = p x ^ } i^o' )■ ^ * s enough to control the first order approximation, 



U) 



exp <^ -Tj[£] 



B^(a^) 



The variables Y i} . . . ,Y n are independent once x is fixed. Thus, in view of our 
target, (15.651) . we need to derive an upper bound of order (1 + 0(e)) for 



l i=i J i=i 



£=0 

t B -i 



E ex P i ce E 1^" f II ^ ( e ^ W ) (5.89) 



£=0 



J'=l 



which holds with Pjf s -probability of order 1 - O(e). 

5.7. Good mesoscopic trajectories. A look at (I5.89D reveals what is to be ex- 
pected from good mesoscopic trajectories. First of all, we may assume that it passes 
through the tube G° N (see (15.13D ) of z*. In particular, x e G° N . Next, by our con- 
struction of the mesoscopic chain P^ B , and in view of (13.20D and (13.211) . the step 
frequencies, Tj[£]/£, are, on average, proportional to pj. Therefore, there exists a 

constant, C\, such that, up to exponentially negligible P^ s -probabilities, 

tj[*b] 



max ■ , r 



(5.90) 



holds. 

A bound on microscopic moment-generating functions. We will now use the 
estimate (15.901) to obtain an upper bound on the product terms in (15.891) . Clearly, 



1 + 0(e)) Mj, uniformly in j and crj j) . Thus, by (T5788T) . 



3 igA,- 



x o(j). 



YjQ- + 0(e)) = — J2(l- e 2 ^) (l {CT(i) =- 1} - tfJPW) = - 1 )) = Y r 



Now, for any t > 0, 



In 



x oU) (ay. 



max V 



a=oO')> s 



2M?"/<r v ^ 



ieAw 



2/i, 



1 



W0=-i} » 



(5.92) 
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where W^°^'' s is the variance with respect to the tilted conditional measure, , 
defined through 

t<>:r(f) - ™ \ • (5.93) 



However, fig 3 ^ s (-) is again a conditional product Bernoulli measure on S^ , i.e., 

^f ,s (-) = <g> iw,) [ • I E *w = Nx °m l • (5 - 94) 



where 



Pile, s) = = ^. (5.95) 

By (I5.90D we need to consider only the case s/Mj < C\. Evidently, there exists 
Si > 0, such that, 

S\ < min min minpAe, s) < max max maxpAe, s) < 1 — Si. (5.96) 

j s<CiMj ieAj j a<C\Mj ieAj 

On the other hand, since x G G° N , there exists <5 2 > 0, such that 

S 2 < mm — — — < max — — — < 1 - d 2 . (5.97) 

j Mj j Mj 

We use the following general covariance bound for product of Bernoulli measures, 
which can be derived from local limit results in a straightforward, albeit painful 
manner. 

Lemma 5.2. Let Si > and S 2 > be fixed. Then, there exists a constant, C = 
C(<5i,5 2 ) < oo, such that, for all conditional Bernoulli product measures on Sm, 
Men, of the form 

M / M \ 

l n ■ |X)& = 2M o) , (5-98) 

i=l V fc=l / 

with pi, . . . ,p M G 1 - Si) and 2M G (-M(l - <5 2 ),M(1 - 5 2 )), and /or ah 
1 < k < I < M, it holds that 

|Cov(l fe= _ 1} ;% = _ 1} )|<^. (5.99) 



Going back to Q5.92D we infer from this that 



UtfJP < exp \ (e>) <T ^ \, (5.100) 



uniformly in I — 0, . . . , In- 

Statistics of mesoscopic trajectories. Q5.89D together with the bound (15. 1001) 
suggests the following notion of goodness of mesoscopic trajectories x: 
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Definition 5.3. We say that a mesoscopic trajectory x = (x^ iA , . . . , x 1b ) is good, 
and write x e T AB , if it passes through G° N , satisfies Q5.90D (and its analog for the 
reversed chain) and, in addition, it satisfies 



iR-l 



£ exp <j 0(e) T ~Jf- \ <Pa,b{^ < 1 + 0(e). (5.101) 

By construction (15.651) automatically holds for any x G T AB . Therefore, our 
target lower bound (15.661) on microscopic capacities will follow from 

Proposition 5.4. Let j AB be the mesoscopic flow constructed in Subsections 15. 31 and 
15.41 and let the set of mesoscopic trajectories T AB be as in Definition ^. 3\ Then Q5.64D 
holds. 

Proof. By (I5.49D we may assume that there exists C > such that, for all x under 
consideration and for all £ = —£a, ■ ■ ■ ,£b — 1, 

0A, B (^,^ + i)<e- C£2/7V . (5.102) 
In view of (15.21) it is enough to check that 

^ x p|°( e )El|-} " X ) 4>AB(^,x e+1 )=0(e), (5.103) 

with P^ s -probabilities of order 1 — o(l). Fix 5 > small and split the sum on the 
left hand side of (I5.103D into two sums corresponding to the terms with £ < iV 1 / 2 5 
and £ > N 1 / 2 ' 5 respectively. Clearly, 

Ef-W. (5.104) 

3=1 3 

uniformly in < £ < iV 1 / 2 5 . On the other hand, from our construction of the 
mesoscopic flow f A B , namely from the choice (15.191) of transition rates inside G° N , 
and from the property (I3.33D of the minimizing curve x(-), it follows that there 
exists a universal (e-independent) constant, K < oo, such that 

P f ^' B (max max ^1 > K ] = o(l). (5.105) 

V 3 £>Ni/2S £pj J 

Therefore, up to P^ B -probabilities of order o(l), the inequality 

0(<-)± T ^<0mV±£ r K*0(4, (S.106) 

j=1 3 j=1 3 

holds uniformly in £ > N l l 2 ~ 5 . A comparison with (I5.102D yields (I5.103D . □ 



The last proposition leads to the inequality (I5.66D . which, together the upper 
bound given in 04.621) . concludes the proof of Theorem 11.31 
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6. Sharp estimates on the mean hitting times 

In this section we conclude the proof of Theorem 11.21 To do this we will use 
Equation (I2.12D with A = S[ml\ and B = S[M], where m* Q is a local minimum 
of Fp^N and M is the set of minima deeper than mj. The denominator on the 
right-hand side of (I2.12D . the capacity, is controlled by Theorem 11.31 What we 
want to prove now is that the equilibrium potential, Ha,b{o'), is close to one in the 
neighborhood of the starting set A, and so small elsewhere that the contributions 
from the sum over a away from the valley containing the set A can be neglected. 
Note that this is not generally true but depends on the choice of sets A and B: the 
condition that all minima m of F^ N such that Fp jN {rn) < F^ N (m* Q ) belong to the 
target set B is crucial. 

In earlier work (see [4]) the standard way to estimate the equilibrium potential 
^a,b(c") was to use the renewal inequality h AjB (a) < cap(Bg) anc ^ bounds on ca- 
pacities. This bound cannot be used here, since the capacities of single points are 
too small. We will therefore use another method to cope with this problem. 

6.1. Mean hitting time and equilibrium potential. Let us start by considering a 
local minimum mj5 of the one-dimensional function Fp t N, and denote by M the set 
of minima m such that Fp^ N (m) < Fp )N (rn* Q ). We then consider the disjoint subsets 
A = S[m* ] and B = S[M]', and write Eq. (I2TT2D as 

erSA ^ ' ' me[-l,l] <jES[m) 

We want to estimate the right-hand side of (I6.1D . This is expected to be of 
order Q/j^jtoq), thus we can readily do away with all contributions where Q^n is 
much smaller. More precisely, we choose 5 > in such a way that, for all N large 
enough, there is no critical point z of Fp )N with F f3tN (z) e [F AJV (mJ), Fp )N {m* Q ) + 5], 
and define 

U 5 = {m : F p>N {m) < F p , N (m£) + 5}. (6.2) 
Denoting by Kg the complement of Us, we obviously have 

Lemma 6.1. 

mAWaA*) < Ne-P NS Qp tN (m* ). (6.3) 

meW| aeS[rn] 

The main problem is to control the equilibrium potential hosier) for configura- 
tions a e S[Us}. To do that, first notice that 

Us=Us(m*) |J Ue(m), (6.4) 

where Us(m) is the connected component of Us containing m (see Fig. 16.11) . Note 
that it can happen that U$(m) = Us(m') for two different minima m, m' e M. 
With this notation we have the following lemma. 

Lemma 6.2. There exists a constant, c > 0, such that, 
(i) for every m e M, 

a£S[U s (m)] 

and 
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U s (m 2 ) 



Figure 4. Decomposition of the magnetization space [—1,1]: is 
represented by dotted lines, while the continuous lines correspond 

U s =Us(m* Q )\J meM U s (m). 



(ii) 



(6.6) 



The treatment of points (i) and (ii) is completely similar, as both rely on a rough 
estimate of the probabilities to leave the starting well before visiting its minimum, 
and it will be discussed in the next section. 

Assuming Lemma 16.21 we can readily conclude the proof of Theorem 11.21 In- 
deed, using (I6.5D together with (I6.3D . we obtain the upper bound 



W,N{°)hA,B{°) < 2^M + 0(Q^K)e- /37Vc ) 



Qj3,N{ml) 



ttN 



2f3a(m* ) 



(6.7) 



where a(mj) is given in (I1.19D . On the other hand, using (16.61) . we get the corre- 
sponding lower bound 



Y Vf3,N(v)h A ,B( a ) > Y Y Vf3,N{v) [1 - (1 - h A)B ((j))} 

Y QwW - o(Q p , N (0^ mc ) 



creS 



> 



Qf3,N(m* ) 



7lN 



2(3a(m* ) 



-,{l + o(l)). 



(6.8) 
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From Equation (I1.12D for Qp,n{™>q) and Equation (I1.31D for ca,p(A, B), we finally 
obtain 

= exp (f3N (F/3 } n(z*) - F AiV (m*))) 



27TiV / /^(l-tanh^/^ + fr)))-! ^ | 
X /3|7i|\/l-»(l-tanh 2 (/?K + / i ))) li + 0lijj ' ^ 

which proves Theorem ll.2l 

6.2. Upper bounds on harmonic functions. We now prove Lemma 16.21 giving 
a detailed proof only for (i), the proof of (ii) being completely analogous. This 
requires, for the first time in this paper, to get an estimate on the minimizer of the 
Dirichlet form, the harmonic function h A ^ B (cr). 

First note that, since h A:B (o~) = P<t(ta < Tb) for all a A U B, the only non zero 
contributions to the sum in (i) come from those sets Us(m) (at most two) whose 
corresponding m is such that there are no minima of M between and m. By 
symmetry we can just analyze one of these two sets, denoted by Us(m*), assuming 
for definiteness that mj < m*. 

Note also that since h AjB (<y) = for all a such that m* < m(a), the problem can 
be reduced further on to the set 

Kg = U 5 (m*) n{m:m< m*}. (6.10) 

Define the mesoscopic counterpart of , namely, for fixed m* G M and n G N, let 
m* G be the minimum of Fp^ N (x) correspondent to m*, and define 

Us = U 5 (m*) = {x G T n N : m(x) G Uj}. (6.11) 

We write the boundary of Us as dUs = d A Us U ObUs, where ObUs = OUs n B, 
and observe that, for all a G S[U$] 

h A)B (a) = ¥ a [r A < t b ] < K[rs[a A u t ] < r S[dBUe] ). (6.12) 
Let max^ pa -C 9(e) <C 1, and for 9 = 9(e) define 

ftS {» effl :E^4). (6.13, 

As before, we denote by OGe the boundary of Go, and write OGe = d A Go U d B Ge, 
where &bGq = dG e n B (see Fig. Ob . 

The strategy to control the equilibrium potential, P ct (ta < t b ), consists in esti- 
mating the probabilities P a [r A < T S[dAGe ] UB }, for tr G «S[C/ 5 \ G 9 ], and P CT [r 5[aAGe] < 
re], for a G G0, in order to apply a renewal argument and to get from these esti- 
mates a bound on the probability of the original event. 

Proceeding on this line, we state the following: 

Proposition 6.3. For any a G (0, 1), there exists n G M, such that the inequality 

K(T A < T S[dAGe]uB ) < e -a-^N[F gM ml)+S-F,Mm(<r))} (6J4) 

holds for all a G S[Us \ Go], n > n , and for all N sufficiently large. 
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Proof of Proposition 16.3b Super-harmonic barrier functions. Throughout the 
next computations, c, d and c" will denote positive constants which are indepen- 
dent on n but may depend on (3 and on the distribution of h. The particular value 
of c and d may change from line to line as the discussion progresses. 
We first observe that, for all a G S[U S \ G e ], 

W*[ta < r s[dA G g ]uB} < ^a[r S [d A u s ] < T S [a A a e ]uB\- (6-15) 

The probability in the r.h.s. of (16.151) is the main object of investigation here. The 
idea which is beyond the proof of bound (16.141) is quite simple. Suppose that ip is 
a bounded super-harmonic function defined on S[U$ \ Ge], i.e. 

(Lip)(a)<0 for all a eS[U s \ G e \. (6.16) 

Then ip(a t ) is a supermartingale, and T = T S[dAUs] A T S[dAGg]uB is an integrable 
stopping time, so that, by Doob's optional stopping theorem, Vct g S[Us \ Ge], 

MW < (6-17) 

On the other hand, 



and hence 



E ff ^(a T ) > mm ip(a')F <T (Ts[d A u s ]<Ts[d A G e ]uB), (6.18) 

K(r s[dAUs] < r s[dAGe]uB ) < ±^- y (6.19) 

The problem is to find a super-harmonic function in order to get a suitable bound 
in (f6TT9l) . 

Proposition 6.4. For any a G (0, 1), there exists n G N such that the function 
= 4>(m(a)), with : IR n i— > R defined as 

(f>(x) = e (i-<*WFpM»), (6.20) 

is super-harmonic in S[Us \ Ge] for all n > n and N sufficiently large. 



METASTABILITY IN THE RFCW MODEL 



46 



(6.21) 



The proof of Proposition 16.41 will involve computations with differences of the 
functions Fp >N . We therefore first collect some elementary properties that we will 
use later. First we need some control on the second derivative of this function. 
From (I3.14D we infer that 

Thus all the potential problems come from the function I N j. 
Lemma 6.5. For any y e (-1, 1), 

tanh _1 (y) - (3e < I' N/ (y) < tanh _1 (y) + /3e, (6.22) 
In particular, as y — > ±1, I' Ne (y) — > ±<x>. 

Proof. Recall that I' N/ (y) = U' N j(y). Set I' N/ (y) = t. Then 

y = jh tanh ^ + P^i) (6.23) 

and hence 

tanh(t -(3e)<y< tanh(t + /?e), (6.24) 

or, equivalently, (16.221) . which proves the lemma. □ 

Lemma 6.6. For any y e (—1, 1) we have that 

- J ^ (y) - i n i-u i^" (6 - 25) 

In particular, for all y e [—1 + z/, 1 — i/], with v e (0, 1/2), 

° £ ^ 5 2, + J + 0(e) £ C ' (6 ' 26) 
and, for all y e (-1, -1 + 1/] U [1 - z/, 1), 

< /^(y) < 7— ft- (6-27) 

1 -\y\ 

Proof. We consider only the case y > 0, the case y < is completely analogous. 
Using the relatio n It e (x) = (U'^ e (I' Ne (x))) 1 and setting te = I' Ne (y)arctanh(y), 
and using Lemma 16.51 we obtain 

^nAv) = — 



EieA^Cl-tanh^^ + t,)) 

< 



|A|(a»)| 

1 



< 



< 



1 - tanh (e/3 + 
1 



1 - tanh 2 (tanh -1 (y) + 2e(3) 
1 

1 - (y + 2e/3tanh / (tanh- 1 (y))) 1 
1 

l-(y + 2 e /?(l-jr>)) 2 ' 



(6.28) 
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where we used that tanh is monotone increasing. The remainder of the proof is 
elementary algebra. □ 

Let us define, for all m such that xg/pe e [—1, 1 — 2/N], 

g e (x) = f (F NtP (x + e t ) - F N ^(x)) . (6.29) 

Lemma 16.61 has the following corollary. 

Corollary 6.7. (i) Ifx e /p e e [—l + u,l — u], with v > 0, then 

gt{x) =-x-h £ + y N/ (x e /pi) + 0(1/ N). (6.30) 

(ii) Ifxt/pt G [-1, -1 + u] U [1 - v, 1 - 2/N], then 

gi(x) = -x - hi + ±I' N j(x e /pi) + 0(1), (6.31) 

where 0(1) is independent of N,n, and v. 
(hi) If xi/pi G [—1 + p, 1 — v\, with v > 0, then there exists c < oo, independent 
of N, such that 

Q 

\gi(x) -gi(x - e £ )\ < — . (6.32) 

(iv) Ifxe/pe G [-1, -1 +v] U [1 - v, 1 - 2/N], then 

\g t (x) - g e (x - e t )\ <C, (6.33) 

where C is a numerical constant independent of N, n, and v. 

The proof of this corollary is elementary and will not be detailed. 
The usefulness of (ii) results from the fact that \V N i/\ is large on that domain. 
More precisely, we have the following lemma. 

Lemma 6.8. There exists v > 0, independent of N and n, such that, ifxg/pi > 1 — v, 
then gt(x) is strictly increasing in x^ and tends to +oo as Xi/p e f +1; similarly if 
xe/pe < —1+u, then ge(x) is strictly decreasing in x e and tends to — oo as xg/pc [ —1. 

Proof Combine (ii) of Corollary 16. 71 with Lemma [6751 and note that h t is bounded 
by hypothesis. □ 

The next step towards the proof of Proposition 16.41 is the following lemma. 

Lemma 6.9. Let m e Ug\Ge and denote by S(m) = {£ : me/pe ^ 1}- Then there 
exists a constant c = c(/3, h) > 0, independent of N and n, such that the following 
holds. If 

£ Pi < g, (6-34) 

then 

e 2 

£ Pi(gi(m)f > c— , (6.35) 

Proof. From the relation I' Ne (x) = U'^\(x), we get that, for all £ e S(m), 

me = Jj Yl tanh V 3 (9e(m)(l + o(l)) + m + hi)) . (6.36) 
Here o(l) tends to zero as N — » oo. 
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We are concerned about small ge(m). Subtracting J2 ieAe tanh (/3 (m + ft*)) on 
both sides of (16.361) and expanding the right-hand side to first order in ge(m), and 
then summing over t e S(m) , we obtain 

1 N ( 1 \ 

- — y^tanh (/3 (m + hi)) - ^ I m e - — ^ tanh (/3(m + ft,;)) J 

i=l e<£S(m) V «eA £ / 

<c P^( m )l< c X] ^( m ) • (6.37) 
ees(m) \£es(m) J 

Notice that the function m i— > m— Xlili tanh (/J (m + ^)) has, by (11.20D . non-zero 
derivative at m* . Moreover, by construction, m* is the only zero of this function in 
Uj(m*). From this observations, together with (I6.37D , we conclude that 

\ 1/2 

^p^ 2 (m)J >c\m-m*\-2 ^ Pt> (6.38) 

for some constant c < oo. Here we used the triangle inequality and the fact that 
\ m t ~ ~k SieAf tann (f3(m + ^t))| < 2p^. Under the hypothesis of the lemma, this 
gives the desired bound if \m — m*\ > d'e/y/6 for some constant c" < oo. On the 
other hand, we can write, for £ e S(m), 

\m e — m}\ < — ^ [tanh ((3 (ge(m)(l + o(l)) + m + fej)) — tanh (f3 (m + 

+ l tann (/^( m + ^)) -tanh(/3(m* + hi))\ 

< cpi\m — m*\ + c' pe\gi(m)\. (6.39) 
Hence we get the bound 

v 1/2 / \ 1/2 



.ees(m) I \ees(m) 



c\m — m 



\ i/ 2 



\ ^ (m £ -m*)^ _ \ ^ 
i. -> Pi i. -> 



Pi 

=1 ££S(m) 



c\m — m 



( \ in 

> c\e 2 /d-A ^ Pi] -c'\rn-m*\ 

\ tyS{m) ) 

> ce/VW -c'\m-m*\ (6.40) 

where in the last line we just used that m ^ G e . The inequalities 06.381) and (16.401) 
now yield (I6.35D . concluding the proof of the lemma. □ 

Proof of Proposition 16. 41 Let o e S[Us \ G e ] and set x = m(a), so that, for ip as in 
Proposition 16.41 Ltjj(a) = L<f)(x). Let a 1 be the configuration obtained from a after 
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a spin-flip at i, and introduce the notation 

n 

L<f ) (x) = Y J L e <p(x), (6.41) 

1=1 

where 

L £ (j)(x) = ^ Pn(v,v 1 )[4>(x + e t ) - <f>(x)] + p^(<T, cr J ) [</>(# - e t ) - <p(x)]. 

(6.42) 

Notice that when x e /p e = ±1, then Af(x) = and the summation over Af(x) in 
A6.42D disappears. 
We define the probabilities 

Y* ±A = PNfao*), (6-43) 

and observe that they are uniformly close to the mesoscopic rates defined in (I4.2D , 
namely 

P CT 

< e C£ , (6.44) 



r N (x,x ± q 
for some c > and e = 1/n. Notice also that 

cpi < P+ i£ + ¥L A < c'p e . (6.45) 

With the above notation and using the convention 0/0 = 0, we get 

L e( f>(x) = <p(x)Fl A [exp (2/5(1 - a)g e {x)) - 1] 

+0(:r)P^ [exp (-2/3(1 - a)^(aj - e t )) - 1] 

= </>(*) (t n ^ (} F° +/ G+(x) + l^^jP^^)) (6.46) 

where we introduced the functions 

G+{x) = exp (2/3(1 - a)^(x)) - 1 + (exp (-2/3(1 - a)^(x - e*)) - 1) (6.47) 

= exp (-2/3(1 - a)gi{x - e t )) - 1 + (exp (2/3(1 - a)^(aj)) - 1) (6.48) 

Ifxi/pi = ±1, the local generator takes the simpler form 

r ^ _ / <K*) F ^ t ex P (-2/3(1 - <*)&(* - e*)) - 1] if ^/p, = 1 

^0(xj - | 0(a . )P ^ [exp (2/3(1 _ _ ^ tf ^ = _j (0.4Vj 

From Lemma [oTBI and inequalities (I6.45D . it follows that, for all £ such that xg/pi = 

±1, 

L t (f>{x) < -(1 + o{l)) P 4{x). (6.50) 
Let us now return to the case when x is not a boundary point. By the detailed 
balance conditions, it holds that 

r N (x, x + e e ) = exp (-2(3g l (x))r N (x + e e , x) , . 

rjv(a:, x - e*) = exp (2(3g e (x - e e ))r N (x - e e , x), 

which implies, together with A6.44D , 

v 

exp (-2f3ge(x) - ce) < =p- < exp (-2fig e (x) + ce) 

-•V (6.52) 
exp (2f3g e (x - e e ) - ce) < ^ < exp (2(3g e (x - e e ) + ce) 
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Inserting the last bounds in (16.471) and (16.481) , and with some computations, we 
obtain 

Gj(x) < (exp (2/3(1 -a)g e {x)) - 1) (1 - exp {2(3ag e {x - e t ) Tee)) (6.53) 
+ exp (2(3gi(x - e e ) =F ce) (exp 2/5(1 - a) (ge(x) - g £ (x - e e )) - 1) 

Gj{x) < (exp (-2/3(1 -a)g e (x-e e ))-l){l- exp (-2/3ag e {x) t ce)) (6.54) 
+ exp (-2/3g e (x) =F ce) (exp 2/3(1 - a) (g e (x) - g e (x - e e )) - 1) 

where =f = -sign (g e (x)) = -sign (g e (x - e e )). 

For all £ such that x e /p e e [—l + u,l — u],we can use (16.321) to get 

Gf(x) < (exp (2/3(1 - a)g E (x)) - 1) (1 - exp {2aPg t {x) =F ce)) + c/iV (6.55) 

G7(cc) < (exp (-2/3(1 - a)g e (x)) - 1) (1 - exp (-2aPg e {x) t ce)) + c/N. (6.56) 
The right hand sides of both (16.551) and (I6.56D are negative if and only if \g e \ > 
Let us define the index sets 

S< = {£: x t /p t e [-1 + - u), \g t (x)\ <§} (6.57) 

S> ee {^laj^et-l + i/,!-!/],!^)^^}. (6.58) 

If £ G S' < , we get immediately that 

max{G+(a;),G7(a;)} < % 2 , (6.59) 

and thus, from (16.46D and (16.451) . 

^0(x) < c -e 2 P 4{x). (6.60) 

To control the r.h.s. of (16.55D and (16.561) when £ e set 

y e ee min {/3 | , ±} < \g e (x) \ . (6.61) 

If g t (x) > then 

exp (2/3(1 - a)^(aj)) - 1 > exp (2(1 - a) w ) - 1 > 2(1 - a)y e (6.62) 

and 

1 - exp (2/3agi(x) - ce) < 1 - exp (ay e ) < -ay e , (6.63) 

so that the product in the r.h.s. of (16.551) is bounded from above by —2(1 — a)ayj. 
On the other hand, if ge(x) < 

exp (2/3(1 - a)g t (x)) - 1 < exp (-2(1 - a)y t ) - 1 < -(1 - a)y t (6.64) 

and 

1 - exp (2(3ag i (x) + ce) > 1 - exp (-aye) > jay e , (6.65) 

and the product in the r.h.s. of (16.551) is bounded from above by — 1(1 — a)ay\. 
Altogether, this proves that, for all £ <E S > , 

Gj(x)<-l(l-a)ayl (6.66) 

and with a similar computation, that 

Gj(x) < -|(1 -a)ay 2 e . (6.67) 

If £ e S > , then we have 

L e 4>(x) < -capey]<f){x). (6.68) 
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It remains to control the case when x e /p e e (—1,-1 + v\ U [1 — u, 1). From 
Lemma |6T8l it follows that, while the positive contribution to Gj(x) and Gj(x) 
remains bounded by a constant, the negative contribution becomes very large as 
soon as v is small enough. More explicitly, for all v small enough, we have 

G$(x) < -(exp(±C") - l) 2 + exp(±C7')(exp(2/?(l - a)c) - 1) < -(1 + o(l)) 
Gj(x) < -(1 - exp(TC')) 2 + exp(TC"){exp{2(3{l - a)c) - 1) < -(1 + o(l)) 

(6.69) 

where C" and C" are positive constants tending to +oo as v \ 0, and the sign ± is 
equal to the sign of x e . Together with (I6.45D and G6.46D , we finally get 



Li<j>{x) < -(1 + o{l))p £ <t>{x). (6.70) 
From (16.501) . (16.601) . (16.68D and (I6.70D . it turns out that the positive contribution 



to the generator L<fi(x) = Y^i=i Li4>(x), comes at most from the indexes £ e S < , 
and can be estimated by 

^ 2 ^. (6.71) 

e&s< 

Now we distinguish two cases according to whether the hypothesis of Lemma 



6.91 are satisfied or not. 

Case 1: J2e^s(x) Pe> i- B Y (16-501), we get 



J2 L ^( x ) < L ^ x ) + L ^ x ) (6 - 72) 

< ~(1 + o(m(x) + ie\ 

which is negative as desired if 9 is small enough, that is, with our choice, if e is 
small enough. 

Case 2: Y^tgs(x) P £ — fe- ^ n tn ^ s case > the assertion of Lemma |6.9| holds. 
By (16301) . (167581) . and (f670l) . we have that, for all £ e S(a;) \ L<, 



Li4>(x) < -p e <p(x) mm{cay 2 , 1} < -captyj(f)(x), (6.73) 
where the last inequality holds for a < 4/c. Let us write the generator as 

i4(x) < Yl L ^ + S (6 - 74) 

The first sum in (I6.74D is bounded from above by 

-ca<j)(x) PfVl ^ -ca(j){x) Y Pe min {P 2 9e(x); \) 

ies(x)\s< ees(x)\s< 

< -ca<P{x) min J /3 2 J] P^ 2 (*); \ \ . (6.75) 
I ees(x)\s< I 



But from Lemma [6791 we know that, for all x eUs\ G$, 



Pe9K^>cj-^s 2 >c"j, (6.76) 
ees(x)\s< 



METASTABILITY IN THE RFCW MODEL 



52 



where c" is a positive constant provided that a > c0. Taking n large enough, it 
holds that 

mini/? 2 PiSti^i) >rmnL"j;lX=c"j, (6.77) 



and then, from (16.71D and (I6.75D , we get 



Lij)((j) < -e 2 {l - a)<p(x)(c aQ- 1 - e'er 1 ). (6.78) 

By our choice of 6 and taking n large enough, the condition d'aO^ 1 — dor 1 > 
a > c6 is satisfied for any a e (0, 1). Hence, for such n's and for N large enough, 
we get that Lij){o) = L<p(x) < concluding the proof of Proposition 16. 4i □ 



Substituting the expression of the super-harmonic function (16.20D in (16.191) , and 
together with (16.151) . we obtain that, for all a 6 S[Us \ G$], 

K[r A < r s[dAGe]uB ] < max e-WWfrM^WMM*))) 

a'&S[d A U s ] 

< e -^-»)l3N[F^ N (m*)+5-F^ N (m(a))]^ (6.79) 

where the last inequality follows from the definition of Us together with the bounds 
in (I3.32D . This concludes the proof of Proposition 16.31 

Renewal estimates on escape probabilities. Let us now come back to the proof 
of Lemma 16.21 An easy consequence of Eq. (I6.14D is that, for all a e S\c)aGq\, 



K(TA < T S[dAGe]uB ) < e ~(^)^MO+S) max e (l-a)PN F8 , N (m) 7 (6 8Q) 

medAGg 

while obviously F a (r A < Ts[d A G e ]uB) = for all a e S[Gg \ OaGo]. To control the 
r.h.s. of fl6.80D . we need the following lemma: 

Lemma 6.10. There exists a constant c < oo, independent of n, such that, for all 
m e Gg, 

Ff3, N (m) < Fp iN (m*) + ce. (6.81) 
Proof. Fix m <E Ge and set m — m* = v. Notice that, from the definition of Gg, 

\\v\\\ < max p e J2 {m£ ~ ^ < e 2 . (6.82) 

Using Taylor's formula, we have 

F PjN {m) = F p>N {m*) + l - {v, A(m*)v) + ^D 3 F^ N (x)v 3 , (6.83) 

where A(m*) is the positive-definite matrix described in Sect. 13.21 (see Eq. (I3.16D ) 
and a; is a suitable element of the ball around m* . From the explicit representation 
of the eigenvalues of A(m*), we see that ||A(ra*)|| < ce^ 1 , and hence 

(v,A(m*)v) < ce^WvWl < ce. (6.84) 



METASTABILITY IN THE RFCW MODEL 



53 



The remainder is given in explicit form as 

D*Fp jN (x)v 3 = E^W^-flE^^^ (6 - 85) 



-y 



1 |A< I" 1 E tgA , tanh(t f + Phj){\ - tanh 2 (t, + gfe)) 3 

3 



^ <=i Pl (|A,|-i £.^(1 - tanh 2 (t, + Ph^y 

where t t = I' Ne (xi/pi). Thus 

n 1 

|D 3 F AJV (x)t; 3 | < cV -jwf < c'e^lHI 2 . < c'e, (6.86) 

where we used that \vt/pt\ < 1. Hence, for some c < 00, independent of n, 

Fp, N {m) < Fp, N {m*) + ce (6.87) 
which proves the lemma. □ 

Inserting the result of Lemma [6.101 into (16.801) . and recalling that Fp^ N (m*) = 

F/3,N(m*), we get that for all a G S[d A Gg] 

P CT (TA < T S[dAGe]uB ) < e -^m(F pM mo)+6-F pM mn-ce)_ (fi gg) 

The last needed ingredient in order to get a suitable estimate on f^ijA < t b ), is 
stated in the following lemma. 

Lemma 6.11. For any 5 2 > 0, there exists n G N, such that, for all n > n , for all 
a G S[dAGg], and for all N large enough, 

K(T B <T S[dAGe] )>e- N ^. (6.89) 

Proof. Fix a 6 S[d^Gg] and set m(0) = ra(cr). As pointed out in the proof of 
Lemma 16.101 every m(0) G d^Gg can be written in the form m(0) = m* + v, 
with v G such that ||t>|| 2 < e. Then, let m = (m(0),m(l), . . . ,m(\\v\\iN) = 
ra*) be a nearest neighbor path in from m(0) to m* , of length A/'|| , u|| 1 , with 
the following property: Denoting by £ t the unique index in {1, . . . ,n} such that 
mg t (t) 7^ rri£ t (t — 1), it holds that 

m £t (t) = m lt (t-l) + |s t , Vt > 1, (6.90) 

where we define 

s t = sign (m| 4 - m £i (t - 1)) . (6.91) 

Note that, by property A6.90D , m(t) G G e for all t > 0. Thus, all microscopic paths, 
(a(t)) t > , such that er(0) = cr and m(a(t)) = m(t), for all t > 1, are contained in 
the event {t b < Ts[d A G e ]}- Thus we get that 
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K(tb < rs[d A G e ]) > K(m( ( r(t)) = m(t),\ft = l,...,\\v\\ 1 N) 

H|iiv 

= ] [ F a (m(a(t)) = m(t)\m(a(t - 1)) = m(t - 1)) 

\\v\\iN 

= n e^w*- 1 )'^- 1 ))- (6 - 92) 

Note that is the set of sites in which a spin-flip corresponds to a step from 

m(t — 1) to m(t). 

The sum of the probabilities in the r.h.s. of (I6.92D corresponds to the quantity 
defined in <H6A3\) . From the inequalities (167441) and (147151) , it follows that, 
for some constant c > depending on (3 and on the distribution of the field, 

> c\A? t (rn(t - 1))\/N > c |A£(m*)|/tf, (6.93) 

where the second inequality follows by our choice of the path m. Now, since 
\Af(m*)\/N = ^(pe±m* e ), using the expression (13.20D for m| and continuing 
from (I6.93D , we obtain 

Inserting the last inequality in (16.92D , and using that, by definition of the path m, 
the number of steps corresponding to a spin-flip in A e is equal to \v e \N, for all 
1 = {l,...,n} , we get 

K{r B < r S [d A Ge]) > n c '^* 



i=l 

n 

= e \\v\\iNWcOTl p \vi\N 

1=1 

> e iV £ ln( C ') e -iV(E? = i^/w) 1/2 ^ 1/2 



> g-^yrv^v, (6.95) 

where in the third line we used the inequality ||v||i < e^ 1 / 2 )!^!^ < -y/e, and in the 
last line we used that ra(0) = m* + u e G$. By our choice of 6 3> e, there exists 
n e N such that, for all n > n , a/| — v^l n (c') < /5<5 2 - For such n's, inequality 
(I6.95D yields the bound (16.891) and concludes the proof of the Lemma. □ 

We finally state the following proposition: 
Proposition 6.12. For all a e S[U$] it holds that 
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Proof. Let us first consider a configuration a G S[d A Gg\. Then it holds 

P<t(ta < tb) < P ct (ta < r 5 [a A £3 e ]us) + P(J ( rA < Ts ' ^ - r "S[^G e ]uAuB) 

< P CT (r A < r 5 [ Sj4Gfl ] UB ) + max ¥ v (r A < r B )¥ a (T S[dAGg] < t b ) 

< K(r A < r s[dAGe]uB ) + max P r ,(r A < r B ) (l - e^ 2 ) , 

ri£S[d A G g ] 

(6.97) 

where in the second line we applied the Markov property, and in the last line we 
insert the result (16.121) . Taking the maximum over a G S[8 a Gq] on both sides of 
(16.971) , and rearranging the summation, we get 

max F a (r A <T B ) < max ¥ a (r A < T S[aA G e ])e pm2 

a&S[d A Gg] a&S[d A G e UB] 

< e -/3JV((l-a)(i^ )2V (m5)+«-F 1 g j ^(m*)-ce)-(52) (6.98) 

where in the last line we used the bound 06.881) . This concludes the proof of (16.961) 
for a G S[d A Gg]. 

Then, let us consider a G S[Ug \ d A G g ]. As before, it holds 

P ct (ta < T B ) < Px(r4 < r s[dAGg]uB ) + ^ P CT (r A < r B , r,, < t s[9aG()]uAuB ) 

v eS[d A G g ] 

< F a (r A < T S[dAGg]uB ) + max P^(r A < T B )F a (T S[dAGg] < t b ) 

ri£S[d A Gg] 

< K{t a < T S[dAGg]uB ) + max F v {t a <t b ), (6.99) 

where P (7 (r j4 < ts^g^ub) is for all a G S[G e \ d A G e ], and exponentially small 
in N for all a e S[U 5 \ G e ] (due to Proposition [6J2). Inserting the bound 06.981) 
in the last equation, provides Eq. (16.96D for a G S[U$ \ d A G g ] and concludes the 
proof. □ 

The proof of formula (I6.5D now follows straightforwardly. From (16.96D , we get 

^ /i A Ar(cr)P <T (r 4 < r s ) 

o-e5[W 4 (m*)] 

< e -PN[(l-a)(F 0>N (m*)+5-F 0>N (m*)-ce)-5 2 ] ^ Q p N (m) 

_ Q p N ( m *ym[aF^ N ( m *)-(l-a)(8-F^ N {m*)-ce)+5 2 } ^ e -0NF PtN (m) 

meUs 

< Q /3iiV (m*)Ar n e /37V [ a ( F ' 3 ' Jv(mS) - F ' 3 ' Jv(m * ) )- (1 ~ a)(,5 " CE)+ ' 52 ], (6.100) 

where in the second inequality we used the expression Q1.9D for Q^^n^o), while 
in the last line we applied the bound Fp iN (m) < Fp :N (m*) = Fp iN (m*), and then 
bounded the cardinality of Us by N n . Finally, choosing a small enough, namely 

a < . ./""T^ ■ (6.101) 

we can easily ensure that (I6.100D implies (I6.5D . 

In exactly the same way one proves (I6.6D . This concludes the proof of Lemma 
Eland thus of Theorem [L2l 
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